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Abstract: This is a review article about neutrino mass models, particularly see-saw mod- 
els involving three active neutrinos which are capable of describing both the atmospheric 



'nI" . neutrino oscillation data, and the large mixing angle MSW solar solution, which is now 

, uniquely specified by recent data. We briefiy review the current experimental status, show 

! how to parametrise and construct the neutrino mixing matrix, and present the leading 

order neutrino Majorana mass matrices. We then introduce the see-saw mechanism, and 
I discuss a natural application of it to current data using the sequential dominance mecha- 

• nism, which we compare to an early proposal for obtaining large mixing angles. We show 

how both the Standard Model and the Minimal Supersymmetric Standard Model may be 
extended to incorporate the see-saw mechanism, and show how the latter case leads to 
, the expectation of lepton fiavour violation. The see-saw mechanism motivates models with 

additional symmetries such as unification and family symmetry models, and we tabulate 
^ I some possible models, before focussing on two particular examples based on 50(10) grand 

^ I unification and either U{1) or SU{3) family symmetry as specific examples. The article 

contains extensive appendices which include techniques for analytically diagonalising dif- 
ferent types of mass matrices involving two large mixing angles and one small mixing angle, 
to leading order in the small mixing angle. 
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1. Introduction 

1.1 Overview for the non-specialist 

In 1930, the Austrian physicist Wolfgang Pauli proposed the existence of particles called 
neutrinos as a "desperate remedy" to account for the missing energy in a type of radioac- 
tivity called beta decay. He deduced that some of the energy must have been taken away 
by a new particle emitted in the decay process, the neutrino. Since then, after decades 
of painstaking experimental and theoretical work, neutrinos have become enshrined as an 
essential part of the accepted quantum description of fundamental particles and forces, the 
Standard Model of Particle Physics. This is a highly successful theory in which elemen- 
tary building blocks of matter are divided into three generations of two kinds of particle - 
quarks and leptons. It also includes three of the fundamental forces of Nature, but does not 
include gravity. The leptons consist of the charged electron, muon and tau, together with 
three electrically neutral particles - the electron neutrino, muon neutrino and tau neutrino. 

Unlike the case for quarks and charged leptons, however, the Standard Model predicts 
that neutrinos have no mass! This might seem curious for a matter particle, but the 
Standard Model predicts that neutrinos always have a left handed spin - rather like rifle 
bullets which spin counter clockwise to the direction of travel. If right-handed neutrinos 
were to be added to the Standard Model, then neutrinos could have the same sort of 
masses as the quarks and charged leptons, and the theory would also predict the existence 
of antineutrinos. However, even without right-handed neutrinos, neutrinos with mass are 
possible, providing that the neutrino is its own antiparticle. Such a mass is then called 
a Majorana mass, named after the Sicilian physicist, Ettore Majorana. But the current 
Standard Model forbids such Majorana masses. These subtle theoretical arguments about 
the nature of neutrinos have now come to the fore, as the results from experiments detecting 
neutrinos from the Sun, as well as atmospheric neutrinos produced by cosmic rays, suggest 
that they do have mass after all. 

The first clues came from an experiment deep underground, carried out by an American 
scientist Raymond Davis Jr., detecting solar neutrinos. It revealed only about one-third 
of the number predicted by theories of how the Sun works. The result puzzled both 
solar and neutrino physicists. However, some Russian researchers, Mikheyev and Smirnov, 
developing ideas proposed previously by Wolfenstein in the U.S., suggested that the solar 
neutrinos might be changing into something else. Only electron neutrinos are emitted 
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by the Sun and they could be converting into muon and tau neutrinos which were not 
being detected on Earth. This effect cahed neutrino oscillations as the types of neutrino 
interconvert over time from one kind to another, was first proposed some time earlier 
by Pontecorvo. The precise mechanism proposed by Mikheyev, Smirnov and Wolfenstein 
involved the resonant enhancement of neutrino oscillations due to matter effects, and is 
known as the MSW effect. 

Most recently, the Sudbury Neutrino Observatory (SNO) in Canada spectacularly 
showed this to be the case. The experiment measured both the flux of the electron neu- 
trinos and the total flux of all three types of neutrinos. The data revealed that physicists' 
theories of the Sun were correct after all. The idea of neutrino oscillations had already 
gained support from the Japanese experiment Super-Kamiokande which in 1998 showed 
that there was a deficit of muon neutrinos reaching Earth when cosmic rays strike the upper 
atmosphere. The results were interpreted as muon neutrinos oscillating into tau neutrinos 
which could not be detected. 

Such neutrino oscillations are analagous to coupled pendulums, where oscillations in 
one pendulum induce oscillations in another pendulum. The coupling strength is defined in 
terms of something called the "mixing angle" . Following the SNO results several research 
groups showed that the electron neutrino must have a mixing angle of about 30 degrees, and 
forms a mass state of 0.007 clcctronvolts or greater (by comparison the electron has a mass 
of about half a megaelectronvolt). The muon and tau neutrinos must have a (maximal) 
mixing angle of about 45 degrees, and form a mass state of about 0.05 electronvolts or 
greater. 

Experimental information on neutrino masses and mixings implies new physics beyond 
the Standard Model, and there has been much activity on the theoretical implications of 
these results. An attractive mechanism for explaining small neutrino masses is the so- 
called see-saw mechanism proposed in 1979 by Murray Gell-Mann, Pierre Ramond and 
Richard Slansky working in the U.S., and independently by Tsutomu Yanagida of Tokyo 
University. The idea is to introduce right-handed neutrinos into the Standard Model which 
are Majorana-type particles with very heavy masses, possibly associated with large mass 
scale at which the three forces of the Standard Model unify. The Heisenberg Uncertainty 
Principle, which allows energy conservation to be violated on small time intervals, then 
allows a left-handed neutrino to convert spontaneously into a heavy right handed neutrino 
for a brief moment before reverting back to being a left-handed neutrino. This results in 
the very small observed Majorana mass for the left-handed neutrino, its smallness being 
associated with the heaviness of the right-handed neutrino, rather like a flea and an elephant 
perched on either end of a see-saw. 

An alternative explanation of small neutrino masses comes from the concept of extra 
dimensions beyond the three that we know of, motivated by theoretical attempts to extend 
the Standard Model to include gravity. The extra dimensions are 'rolled up' on a very small 
scale so that they are not normally observable. It has been suggested that right-handed 
neutrinos (but not the rest of the Standard Model particles) experience one or more of 
these extra dimensions. The right handed neutrinos then only spend part of their time in 
our world, leading to apparently small neutrino masses. 
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Cosmology today presents two major puzzles: why there is an excess of matter over 
antimatter in the Universe; and what is the major matter constituent of the Universe? 
Massive neutrinos may hold important clues. 

Matter and antimatter would have been created in equal amounts in the Big Bang 
but all we see is a small amount of excess matter. The see-saw mechanism allows for 
a novel resolution to this puzzle. The idea, due to Masataka Fukugita and Tsutomu 
Yanagida of Tokyo University, is that when the Universe was very hot, just after the Big 
Bang, the heavy right-handed neutrinos would have been produced, and could have decayed 
preferentially into leptons rather than antileptons. The see-saw mechanism therefore opens 
up the possibility of generating the baryon asymmetry of the universe via "leptogenesis" . 
This process requires CP violation for neutrinos which could be studied experimentally by 
firing a very intense neutrino beam right through the Earth and detecting it with a huge 
neutrino detector when it emerges. 

Studies of the movements of galaxies and galaxy clusters suggest that at least 90 per 
cent of the mass of the Universe is made of unknown dark matter. Cosmology is sensitive 
to the absolute values of neutrino masses, in the form of relic hot dark matter. Neutrinos 
could constitute anything from 0.1 to 2 per cent of the mass of the Universe, corresponding 
to the heaviest neutrino being in the mass range 0.05 to about 0.23 electronvolt. Neutrinos 
any heavier than this would lead to galaxies being less clumped than actually observed by 
the recent 2dF Galaxy Redshift Survey. This illustrates the breathtaking rate at which 
neutrino physics continues to advance. 

1.2 About this review 

There are many good reviews already in the literature, for example |^], 0]. Three 

possible ways to extend the Standard Model in order to account for the neutrino mass 
spectrum are the see-saw mechanism extra dimensions [^], and R-parity violating su- 
persymmetry In this review we focus on theoretical approaches to understanding 
neutrino masses and mixings in the framework of the see-saw mechanism, assuming three 
active neutrinos. The goal of such models is to account for two large mixing angles, one 
small mixing angle, and a pattern of neutrino masses consistent with observation. We are 
now in the unique position in the history of neutrino physics of knowing not only that 
neutrino mass is real, and hence the Standard Model at least in its minimal formulation is 
incomplete, but also we have a unique solution to the solar neutrino problem in the form 
of the large mixing angle solution. In this sense a review of neutrino mass models is very 
timely, since it has only been within the last year that the solar solution has been uniquely 
specified. That combined with the atmospheric oscillation data severely constrains theoret- 
ical models, and in fact rules out many possibilities which predicted other solar solutions. 
Of course many possibilities remain, and we shall mention several of them here. However 
this review is not supposed to be an encyclopaedic review of all possible models, but instead 
a review of useful approaches and techniques that may be applied to constructing different 
classes of models. 

We give a strong emphasis to classes of models where the two large mixing angles can 
arise naturally and consistently with a neutrino mass hierarchy, and although we classify 



- 4 - 



all possible neutrino mass structures, we do not spend much time on those structures 
which apparently require a high degree of fine-tuning to achieve. We show that if one 
of the right-handed neutrinos contributes dominantly in the see-saw mechanism to the 
heaviest neutrino mass, and a second right-handed neutrino contributes dominantly to 
the second heaviest neutrino mass, then large atmospheric and solar mixing angles may 
be interpreted as simple ratios of Yukawa couplings. We refer to this natural mechanism 
as sequential dominance. Although sequential dominance looks very specialised it is not: 
either the right-handed neutrinos contribute equally via the see-saw mechanism to neutrino 
masses, or some of them contribute more than others. The second possibility corresponds 
to sequential dominance, and allows a very natural and intuitively appealing explanation 
of the neutrino mass hierarchy with two large mixing angles. Sequential dominance is not 
a model, it is a mechanism in search of a model. The conditions for sequential dominance, 
such as ratios of Yukawa couplings being of order unity for large mixing angles, and the 
required pattern of right-handed neutrino masses are put in by hand and require further 
theoretical input. This motivates models with extra symmetry such as unified models 
and models with family symmetry, which we briefly review. There are a huge number of 
proposals in the literature, but assuming sequential dominance, and the important clues 
provided by quark masses and mixing angles, severely constrains the possible successful 
models. Wc discuss one particularly successful model as an example, but of course there 
may be others, but maybe not so many as may be thought at first. 

The layout of the remainder of the review article is as follows. In section 2 we introduce 
and review the current status of neutrino masses and mixing angles. We also parametrise 
the neutrino mixing matrix in two different ways, whose equivalence is discussed in an 
Appendix. We show how it may be constructed theoretically from the underlying mass 
matrices, and then show how the proceedure may be driven the other way to derive the form 
of the neutrino mass matrix whose leading order forms may be classified. The properties of 
the matix corresponding to hierarchical neutrino masses are explored. Section 3 introduces 
the see-saw mechanism, which is central to this review, in both its simplest version, and 
including more complicated versions. In section 4 we show how the see-saw mechanism may 
be applied to the hierarchical case in a very natural way using sequential dominance, discuss 
different types of sequential dominance, and a link with leptogenesis. We also discuss an 
alternative early approach to obtaining large mixing angles from the see-saw mechanism, 
and show that it is quite different from sequential dominance. Section 5 incorporates the 
see-saw mechanism into the Standard Model, and its Supersymmetric version, where it 
leads to lepton fiavour violation. In section 6 we go beyond these minimal extensions of 
the Standard Model, or its supersymmetric version, and show how the see-saw mechanism 
motivates ideas of Unification and Family Symmetry, and briefiy review the huge literature 
that has grown up around such approaches, before focussing on two particular models 
based on 50(10) grand unification and either U{1) or SU{3) family symmetry. Section 7 
concludes the review. 

We also present extensive appendices which deal with more technical issues, but which 
may provide useful model building tools. Appendix A proves the equivalence between 
different parametrisations of the neutrino mixing matrix and gives a useful dictionary. 
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Appendix B gives the full three family neutrino oscillation formula (in vacuum). Appendix 
C derives the formula given in the text for charged lepton contributions to the neutrino 
mixing matrix. Finally Appendix D discusses in detail how to diagonalise different kinds of 
mass matrices involving two large mixing angles analytically to leading order in the small 
mixing angle. 



2. Neutrino Masses and Mixing Angles 

The history of neutrino oscillations dates back to the work of Pontecorvo who in 1957 Q 
proposed v ^ D oscillations in analogy with K ^ K oscillations, described as the mixing of 
two Majorana neutrinos. Pontecorvo was the first to realise that what we call the "electron 
neutrino" for example is really a linear combination of mass eigenstate neutrinos, and that 
this feature could lead to neutrino oscillations of the kind f e — > f/^ Q. Later on MSW 
proposed that such neutrino oscillations could be resonantly enhanced in the Sun . The 
present section introduces the basic formalism of neutrino masses and mixing angles, gives 
an up-to-date summary of the current experimental status of this fast moving field, and 
discusses future experimental prospects. Later in this section we also discuss some more 
theoretical aspects such as charged lepton contributions to neutrino mixing angles, and the 
neutrino mass matrix. 

2.1 Two state atmospheric neutrino mixing 



In 1998 the Super-Kamiokande experiment published a paper |11] which represents a wa- 
tershed in the history of neutrino physics. Super-Kamiokande measured the number of 
electron and muon neutrinos that arrive at the Earth's surface as a result of cosmic ray 
interactions in the upper atmosphere, which are referred to as "atmospheric neutrinos". 
While the number and and angular distribution of electron neutrinos is as expected, Super- 
Kamiokande showed that the number of muon neutrinos is significantly smaller than ex- 
pected and that the fiux of muon neutrinos exhibits a strong dependence on the zenith 
angle. These observations gave compelling evidence that muon neutrinos undergo flavour 
oscillations and this in turn implies that at least one neutrino flavour has a non-zero mass. 
The standard interpretation is that muon neutrinos are oscillating into tau neutrinos. 

Current atmospheric neutrino oscillation data are well described by simple two-state 
mixing 

( 1 _ [ ^23 sin 6*23 

j y - sin 6*23 cos 6*23 

and the two-state Probability oscillation formula 




(2.1) 



P{v^ Vr) = sin^ 26*23 s\n^{l.27/^ml2L/ E) (2.2) 



where 

Am^j = mf — m| (2-3) 

and rrii are the physical neutrino mass eigenvalues associated with the mass eigenstates fj. 
Am|2 is in units of eV^, the baseline L is in km and the beam energy E is in GeV. 
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Figure 1: Summary of the currently allowed regions from a global analysis of atmospheric and 
solar neutrino experiments including first results from KamLAND (from H.Murayama's web site 
http: / /hitoshi. berkeley.edu / neutrino / . ) 

The atmospheric data is statistically dominated by the Super-Kamiokande results and 
the latest reported data sample as of the time of writing leads to: 

sin^ 26*23 > 0.92, 1.3 x IQ-^eV^ < {Aml^l < 3.0 x W'^eV'^iMCL) (2.4) 

The Super-Kamiokande region is shown in Fig.|l|. The atmospheric neutrino data is thus 
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consistent with maximal — Vr neutrino mixing 622, ~ 7r/4 with |Am|2| ~ 2.5 x 10~'^eV^ 
and the sign of Amig undetermined. The maximal mixing angle means that we identify 
the heavy atmospheric neutrino of mass ms as being approximately 

.3 - (2.5) 

and in addition there is a lighter orthogonal combination of mass m2, 

- (2.6) 

2.2 Three family solar neutrino mixing 

Super-Kamiokande is also sensitive to the electron neutrinos arriving from the Sun, the "so- 
lar neutrinos" , and has independently confirmed the reported deficit of such solar neutrinos 
long reported by other experiments. For example Davis's Homestake Chlorine experiment 
which began data taking in 1970 consists of 615 tons of tetrachloroethylene, and uses ra- 
diochemical techniques to determine the Ar37 production rate. More recently the SAGE 
and Gallex experiments contain large amounts of Ga71 which is converted to Ge71 by low 
energy electron neutrinos arising from the dominant pp reaction in the Sun. The combined 
data from these and other experiments implies an energy dependent suppression of solar 
neutrinos which can be interpreted as due to flavour oscillations. Taken together with 
the atmospheric data, this requires that a second neutrino flavour has a non-zero mass. 
The standard interpretation is that the electron neutrinos oscillate into the light linear 

combination V2 ~ '^'^ /"^ • 
V 2 

SNO measurements of charged current (CC) reaction on deuterium is sensitive exclu- 
sively to i^e's, while the elastic scattering (ES) off electrons also has a small sensitivity to 
z/^'s and i^^'s. The CC ratio is significantly smaller than the ES ratio. This immediately 
disfavours oscillations of v'^s to sterile neutrinos which would lead to a diminished flux of 
electron neutrinos, but equal CC and ES ratios. On the other hand the different ratios 
are consistent with oscillations of v^s to active neutrinos f^'s and z/t-'s since this would 
lead to a larger ES rate since this has a neutral current component. The SNO analysis 
is nicely consistent with both the hypothesis that electron neutrinos from the Sun oscil- 
late into other active flavours, and with the Standard Solar Model prediction. The latest 
results from SNO including the data taken with salt inserted into the detector to boost 



the efficiency of detecting the neutral current events |12], strongly favour the large mixing 
angle (LMA) MSW solution. In other words there is no longer any solar neutrino problem: 
we have instead solar neutrino mass! 

The minimal neutrino sector required to account for the atmospheric and solar neutrino 
oscillation data thus consists of three light physical neutrinos with left-handed flavour 
eigenstates, fe, t'^, and defined to be those states that share the same electroweak 
doublet as the left-handed charged lepton mass eigenstates. Within the framework of 
three-neutrino oscillations, the neutrino flavor eigenstates u^, Vfi, and v-r are related to the 
neutrino mass eigenstates 1^1, 1^2, and with mass mi, m2, and ma, respectively, by a 3 x 3 
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unitary matrix called the lepton mixing matrix ?7 |13] 



/ Uel Ue2 Ue3 

U^i U^2 U^s 

\Url Ur2 Ur3 




(2.7) 



Assuming the light neutrinos are Majorana, U can be parameterized in terms of three 
mixing angles 6ij and three complex phases 5ij. A unitary matrix has six phases but three 
of them are removed by the phase symmetry of the charged lepton Dirac masses. Since 
the neutrino masses are Majorana there is no additional phase symmetry associated with 
them, unlike the case of quark mixing where a further two phases may be removed. 

If we suppose to begin with that the phases are zero, then the lepton mixing matrix 
may be parametrised by a product of three Euler rotations. 



U — R23R13R12 



where 



23 



R 



13 



R 



12 




(2.9) 



(2.10) 



(2.11) 



cos 6ij and Sij 



sin 9ij . Note that the allowed range of the angles is < 9ij < 



where Cj 
7r/2. 

CHOOZ is a reactor experiment that failed to see any signal of neutrino oscillations 
over the Super-Kamiokande mass range. CHOOZ data from ~^ disappearance not 
being observed provides a significant constraint on ^13 over the Super-Kamiokande (SK) 



prefered range of Am|2 



sin^ 26*13 < 0.1 - 0.3 



(2.12) 



The CHOOZ experiment therefore limits Q\3 ^ 0.2 over the favoured atmospheric range, 
as shown in Fig.||. 

KamLAND is a more powerful reactor experiment that measures z/g's produced by 
surrounding nuclear reactors. KamLAND has already seen a signal of neutrino oscillations 
over the LMA MSW mass range, and has recently confirmed the LMA MSW region "in the 
laboratory" |]l5|. KamLAND and SNO results when combined with other solar neutrino 
data especially that of Super-Kamiokande uniquely specify the large mixing angle (LMA) 
MSW fTO] solar solution with three active light neutrino states, a large solar angle 



tan^ 6*12 ~ 0.4, Am^^ « 7 x lO^^eV^. 



(2.13) 
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according to the most recent global fits performed after the SNO salt data |T^] . Kam- 
LAND has thus not only confirmed solar neutrino oscillations, but have also uniquely 
specified the large mixing angle (LMA) solar solution, heralding a new era of precision 
neutrino physics. 

The currently regions of atmospheric and solar parameter space allowed by all experi- 
ments are depicted in Figure |^. ^ In Figure [l| the atmospheric and LMA MSW solar regions 
are clearly shown as elliptical regions, with the SMA, LOW and VAC regions now having 
disappeared. One of the KamLAND rate plus shape allowed regions shown in Figure || 
intersects the central part of the LMA ellipse near the best fit LMA point as determined 
from the solar data alone, thereby confirming the LMA MSM solution. 



2.3 Summary of neutrino mixing angles and mass patterns 

The current experimental situation is summarized by 623 ~ 7r/4, ^13 < 0.2, and 612 « tt/6. 
Ignoring phases, the relation between the neutrino flavor eigenstates u^, i^f^, and z/,- and the 
neutrino mass eigenstates vi, U2, and 1^3 is just given as a product of three Euler rotations 
in Eq.2.8 as depicted in Fig.0. This corresponds to the approximate form of mixing matrix 




Figure 2: The relation between the neutrino flavor eigenstates Ve, v^, and Vt and the neutrino 
mass eigenstates vi, V2, and v^i, in terms of the three mixing angles 9i2, ^13, ^23- The atmospheric 
angle is ^23 ~ 7''/4, the CHOOZ angle is 6'i3 < 0.2, and the solar angle is 9i2 ~ 7r/6. 
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where 9i2 ~ 7r/6 corresponds to S12 
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(2.14) 



It is clear that neutrino oscillations, which only depend on Am; 



information about the absolute value of the neutrino mass squared eigenvalues m^, and 



^For more detailed most up to date plots of the LMA MSW region see 
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Figure 3: Alternative neutrino mass patterns that are consistent with neutrino oscillation expla- 
nations of the atmospheric and solar data. The absolute scale of neutrino masses is not fixed by 
oscillation data and the lightest neutrino mass may vary from 0.0-0.23 eV. 

there are basically two patterns of neutrino mass squared orderings consistent with the 
atmospheric and solar data as shown in Fig.H. 



2.4 Three family neutrino mixing with phases 

Including the phases, assuming the Hght neutrinos are Majorana, U can be parameterized 
in terms of three mixing angles 6ij, a Dirac phase 6, together with two Majorana phases 
Pi, (32, as follows 

U = R23U13RUP12 (2.15) 



where 



U 



13 



12 



/ ci3 Osige-*^' 
1 

\ -5136^^ CIS 

/ e^^i \ 
e^^2 
1 



(2.16) 



(2.17) 



where Cij = cos 6ij and Sij = sin 6ij, and R23, R12 were defined below Eq.|2.8|. 

Alternatively the lepton mixing matrix may be expressed as a product of three complex 
Euler rotations, 

U = U23U13U12 (2.18) 



where 





(1 





U23 = 





C23 






-S23e^'523 





5236 
C23 



(2.19) 
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(2.20) 



(2.21) 



The equivalence of different parametrisations of the lepton mixing matrix, and the 
relation between them is discussed in Appendix A. 

Three family oscillation probabilities depend upon the time-of-flight (and hence the 
baseline L), the Amfj, and U (and hence ^i2)^23j^13) and 6). Three-state neutrino mixing 
is discussed in Appendix B. Since we have assumed that the neutrinos are Majorana, there 
are two extra phases, but only one combination S = (5i3 — 623 — S12 affects oscillations. If 
the neutrinos are Dirac, then the phases /3i = /32 = 0, but the phase 6 remains. 



2.5 The LSND signal 

The signal of another independent mass splitting from the LSND accelerator experiment 



17] would either require a further light neutrino state with no weak interactions (a so-called 



"sterile neutrino") or some other non-standard physics. This effect has not been confirmed 
by a similar experiment KARMEN [18|, and currently a decisive experiment MiniBooNE 
is underway to decide the issue. In Figure |l] the LSND signal region is indicated, together 
with the KARMEN excluded region. 



2.6 Future experimental prospects 

Further experimental progress from SNO and KamLAND will consist of pinning down 
LMA MSW parameters to high accuracy. Neutrino physics has, now entered the precision 
era. Future neutrino oscillation experiments, will give accurate information about the 
mass squared splittings Am?^- = mf — mj, mixing angles, and CP violating phase. In 
the near future much better solar neutrino measurements will be available as KamLAND, 
SNO and Borexino furnish us with new and better data. The K2K long baseline (LBL) 
experiment from KEK to Super-Kamiokande has recently reported results in its phases I 
and II, which cover the atmospheric region and support the Super-Kamiokande results. 
In the longer term LBL experiments such as MINOS and eventually the CERN to Gran 
Sasso experiments will given more accurate determinations of the atmospheric parameters, 
eventually to 10%. J-PARC will be an "off-axis superbeam" over a LBL of 295 km to 
Super-Kamiokande due to start in 2008. Its first goal is to measure ^13 or set a limit on 
it of about 0.05 (as compared to the CHOOZ limit on ^13 of about 0.2). Interestingly 
MINOS over a LBL of 735 km is more sensitive than J-PARC to matter effects, so there 
should be some interesting complementarity between these two experiments, which could 
for example allow the sign of Am^g to be determined. The ultimate goal of oscillation 
experiments however is to measure the CP violating phase 6. An upgraded J-PARC with a 
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4MW proton driver and a 1 megaton Hyper-Kamiokande detector, or some sort of Neutrino 
Factory based on muon storage rings would seem to be required for this purpose [p!g| ]. 

Oscillation experiments are not capable of telling us anything about the absolute scale 
of neutrino masses. The Tritium beta decay experiment KATRIN will tell us about the 
absolute scale of neutrino mass down to about 0.35 eV. The neutrinoless double beta decay 
experiment GENIUS will probe the Majorana nature of the electron neutrino down to 
about 0.01 eV [20|. Recent results from the 2dF galaxy redshift survey and WMAP, when 



combined with oscillation data, give the strong limit on the absolute mass of each neutrino 



species of about 0.23 eV ||2l|, |22|]. Turning to astrophysics, a galactic supernova could give 
valuable information about neutrino masses [^] . In future detection of energetic neutrinos 
from gamma ray bursts (GRBs), by neutrino telescopes such as ANTARES or ICECUBE 
could also provide important astrophysical information, and may provide another means 



of probing neutrino mass, and even quantum gravity [24]. 



2.7 Charged lepton contributions to neutrino masses and mixing angles 

Although we refer to neutrino masses and mixing angles, it is worth pointing out that in 
general they could originate, at least in part, from the charged lepton sector. The (low 
energy) Lagrangian involving the charged lepton and neutrino mass matrices 

(^mass = - {eLieL2eL3)'mf^{eRieR2eR3)'^ 

- ]^{l'LlVL2VLz)mLL{T^liVl2vlr^Y + H.c. (2.22) 

where en are the three left-handed charged lepton states, are the right-handed charged 
lepton states, vii are the three left-handed neutrino states, and ^'£j are their CP conjugates. 
Note that the states un are not the mass eigenstate neutrinos since rriLL is not diagonal 
in general. We shall refer to the mass eigenstate neutrinos as Vi (without the L subscript). 



as in Eq. 2.7 



In general the neutrino and charged lepton masses are given by the eigenvalues of a 
complex charged lepton mass matrix mf^^ and a complex symmetric neutrino Majorana 
matrix niLL-, obtained by diagonalising these mass matrices, 

/ 



I (2.23) 
\ ^ 

/ mi 

m2 I (2.24) 



Y m3 

where V^^ , V^'^ , V^'^ are unitary tranformations on the left-handed charged lepton fields 
El, right-handed charged lepton fields Er, and left-handed neutrino fields vl which put 
the mass matrices into diagonal form with real eigenvalues. 

After having diagonalised the mass matrices, the lepton mixing matrix is then con- 
structed by 

U = V'^W^'^ (2.25) 
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A unitary three dimensional matrix has six independent phases. As discussed in Appendix 
B, the freedom in the charged lepton phase enables three of the phases to be removed from 
U to leave three phases. Since we have assumed that the neutrinos are Majorana, there 
is no further phase freedom, and the three remaining phases are physical (unlike the case 
of Dirac neutrinos where a further two phases can be removed, analagous to the case of 
the CKM matrix in the quark sector). Having constructed the lepton mixing matrix as 
discussed above, it may then be parametrised as discussed in section 2. Having done this 
one may then ask how much of a contribution to a particular mixing angle or phase, comes 
from the neutrino sector, and how much comes from the charged lepton sector. The lepton 



mixing matrix is constructed in Eq. |2.17| as a product of a unitary matrix from the charged 
lepton sector V^^ and a unitary matrix from the neutrino sector y^^. Each of these 
unitary matrices may be parametrised by its own mixing angles and phases analagous to 
the lepton mixing matrix parameters. As shown in Appendix C [^] the lepton mixing 
matrix can be expanded in terms of neutrino and charged lepton mixing angles and phases 
to leading order in the charged lepton mixing angles which are assumed to be small, 

3236-'^^^ ~ S23e-'^"^^ - ^fg^c^f e-^'^S'' (2.26) 

(2.27) 



'12 *23^ 



'12 ^23^12 



as "-12 •'23' 

^-c^fc^'^e-^" (2.28) 



Clearly ^13 receives important contributions not just from ^^'f j but also from the charged 
lepton angles 1 ^13 ■ models where ^J'g is extremely small, ^13 may originate 
almost entirely from the charged lepton sector. Charged lepton contributions could also 
be important in models where ^^'2 — '^/^^ since charged lepton mixing angles may allow 
consistency with the LMA MSW solution. Such effects are important for the inverted 
hierarchy model ||2^ . 

Note that it is useful and possible to be able to diagonalise the mass matrices analyti- 
cally, at least to first order in the small 13 mixing angles, but allowing the 23 and 12 angles 
to be large, while retaining all the phases. The proceedure for doing this is discussed for a 
hierarchical general mass matrix in Appendix D.l, for a hierarchical neutrino mass matrix 
in Appendix D.2, and for an inverted hierarchical neutrino mass matrix in Appendix D.3. 
The analytic results in these Appendices enable the separate mixing angles and phases 
associated with each of the unitary transformations V^'^ and ^ to be obtained in many 
useful cases of interest. 

2.8 The neutrino mass matrix 

For many (but not all) purposes it is convenient to forget about the division between 
charged lepton and neutrino mixing angles and work in a basis where the charged lepton 
mass matrix is diagonal. Then the lepton mixing angles and phases simply correspond to 
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the neutrino ones. In this special basis the mass matrix is given from Eg. 2.24 and Eg, 2. 17 
as 



rriLL = U 



/ mi 
m2 
y ms 



(2.29) 



For a given assumed form of U and set of neutrino masses rrii one may use Eg.2.2£ to 
"derive" the form of the neutrino mass matrix rriLL, and this results in the candidate mass 
matrices in Table |l] ||28[] . Only the leading order forms are displayed explicitly in Table |l|, 
and more accurate structures may be obtained case by case. 



Table 1: Leading order low energy neutrino Majorana mass matrices rriLL consistent with large 
atmospheric and solar mixing angles, classified according to the rate of ncutrinoless double beta 
decay and the pattern of neutrino masses. 





Type I 
Small /3/3oiy 


Type II 
Large PPoi^ 


A 

Normal hierarchy 
ml, m| <C m| 




3^ < O.OC 

/o 0^ 
1 1 


82 eV 

m 
2 




B 

Inverted hierarchy 
ml 3> m§ 




3^ < O.OC 

^01l\ 
1 

^looj 


82 eV 

m 
^/2 




)^ > 0.00 

'^l 0^ 

i i 
"22 


85 eV 

TO 


C 

Approximate degeneracy 

9 ? ? 
TOJ « 7712 ~ '^S 


( ' 717l\ 

J- I I TO 

V2 2 2 

\ J_ 1 1 / 

\ V2 2 2 / 


PP 
c 


0^ ^ O.OJ 
liag(l,l,] 
^1 0^ 
1 


>5 eV 

)m 

TO 



In Table |l| the mass matrices are classified into two types: 
Type I - small neutrinoless double beta decay 
Type II - large neutrinoless double beta decay 
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They are also classified into the limiting cases consistent with the mass squared order- 
ings in Fig.^: 

A - Normal hierarchy nif, <^ m| 

B - Inverted hierarchy ~ ^ m'^ 

C - Approximate degeneracy mf ~ m| ~ m| 

Thus according to our classification there is only one neutrino mass matrix consistent 
with the normal neutrino mass hierarchy which we call Type lA, corresponding to the 
leading order neutrino masses of the form rui = (0,0, m). For the inverted hierarchy there 
are two cases. Type IB corresponding to nii = (m, —m, 0) or Type IIB corresponding to 
nii = (m,m,0). For the approximate degeneracy cases there are three cases, Type IC 
correponding to m-j = (m, — m, m) and two examples of Type IIC corresponding to either 
rrii = (m, m, m) or rui = (m, m, —m). 

At present experiment allows any of the matrices in Table |l|. In future it will be 
possible to uniquely specify the neutrino matrix in the following way: 

1. Neutrinoless double beta effectively measures the 11 element of the mass matrix 
rriLL corresponding to 

(3Po. = Yl ^er^i (2.30) 
i 

and is clearly capable of resolving Type I from Type II cases according to the bounds given 
in Table |l] [^9|. There has been a recent claim of a signal in neutrinoless double beta decay 
correponding to /?/3oi/ = 0.11 — 0.56 eV at 95% C.L. |3C]. However this claim has been 



criticised by two groups |31], [^] and in turn this criticism has been refuted |33]. Since 
the Heidelberg-Moscow experiment has almost reached its full sensitivity, we may have to 
wait for a next generation experiment such as GENIUS |2C] which is capable of pushing 



down the sensitivity to 0.01 eV to resolve this question. 

2. A neutrino factory will measure the sign of ^"^^ resolve A from B. 

3. Tritium beta decay experiments are sensitive to C since they measure the "electron 
neutrino mass" defined by 

:^|C7ei|Vi|- (2.31) 



For example the KATRIN |34| experiment has a proposed sensitivity of 0.35 eV. As already 
mentioned the galaxy power spectrum combined with solar and atmospheric oscillation 
data already limits each neutrino mass to be less than about 0.23 eV, and this limit is also 
expected to improve in the future. Also it is worth mentioning that in future it may be 
possible to measure neutrino masses from gamma ray bursts using time of flight techniques 
in principle down to 0.001 eV |24]. 

Type IIB and C involve small fractional mass splittings |Am^j| <C rin? which are un- 



stable under radiative corrections |35|, and even the most natural Type IC case is difficult 
to implement |3^,||3^. Types lA and IB seem to be the most natural cases. 

Consider the case of full neutrino mass hierarchy ms ^ m2 ^ mi ~ 0, which is 
a special case of Type lA, where in this case ~ Y^|Am|^ ~ 5.10^^ eV and m2 ~ 
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a4 



7.10 eV. From Eqs.2.14,2.29 we find the symmetric mass matrix, 



rriLL 



. i(m3+n.2C?,) 

V • 



-^{m2Si2Ci2 - 1713013) 



7712 cf 2) 



(2.32) 



l(m3 + m2C?2) 



neglecting terms like 17126 13. Clearly this expression reduces to the leading Type I A form 
with m = 1113 in the approximation that 771-2 and ^13 are neglected. However the more exact 



expression in Eq.2.32 shows that the required form of ttt-ll should have a very definite 



detailed structure, which goes beyond the leading approximation in Table |^. For example 
the requirement 1112 <^ 771-3 implies that the sub-determinant of the mass matrix miL is 
small: 

in22 "723 
"123 "^33 



det 



<C 7773. 



(2.33) 



This requirement in Eq. p.33 is satisfied by Eq. p.32| , as may be readily seen, and this condi- 
tion must be reproduced in a natural way (without fine-tuning) by any successful theory. 



3. The See- Saw Mechanism 



There are several different kinds of see-saw mechanism in the literature. In this review 
we shall focus on the simplest Type I see-saw mechanism, which we shall introduce below. 
However for completeness we shall also discuss the type H see-saw mechanism and the 
double see-saw mechanism. 



3.1 Type I See-Saw 

Before discussing the see-saw mechanism it is worth first reviewing the different types of 
neutrino mass that are possible. So far we have been assuming that neutrino masses are 
Majorana masses of the form 

ihllW^^l (3-1) 

where is a left-handed neutrino field and is the CP conjugate of a left-handed neutrino 
field, in other words a right-handed antineutrino field. Such Majorana masses are possible 
to since both the neutrino and the antineutrino are electrically neutral and so Majorana 
masses are not forbidden by electric charge conservation. For this reason a Majorana mass 
for the electron would be strictly forbidden. However such Majorana neutrino masses vio- 
late lepton number conservation, and in the standard model, assuming only Higgs doublets 
are present, are forbidden at the renormalisable level by gauge invariance. The idea of the 
simplest version of the see-saw mechanism is to assume that such terms are zero to begin 
with, but are generated effectively, after right-handed neutrinos are introduced |5|. 

If we introduce right-handed neutrino fields then there are two sorts of additional 
neutrino mass terms that are possible. There are additional Majorana masses of the form 



(3.2) 



where i^r is a right-handed neutrino field and u'^ is the CP conjugate of a right-handed 
neutrino field, in other words a left-handed antineutrino field. In addition there are Dirac 
masses of the form 

rriLBpLi'R- (3.3) 

Such Dirac mass terms conserve lepton number, and are not forbidden by electric charge 
conservation even for the charged leptons and quarks. 



Once this is done then the types of neutrino mass discussed in Eqs.3.2,p^ (but not 



Eq j3.1| since we assume no Higgs triplets) are permitted, and we have the mass matrix 

(3.4) 



Mrr 

Since the right-handed neutrinos are electroweak singlets the Majorana masses of the right- 
handed neutrinos Mrr may be orders of magnitude larger than the electroweak scale. In 



the approximation that Mrr ^ rniR the matrix in Eq.3.4 may be diagonalised to yield 



eff'ective Majorana masses of the type in Eq.3.1 



rriLL = -rriLRMj^l^ni^^. (3.5) 

The eff'ective left-handed Majorana masses rriLL are naturally suppressed by the heavy 
scale Mrr. In a one family example if we take rriLR = M\y and Mrr = Mqut then we 
find rriLL ~ 10~^ eV which looks good for solar neutrinos. Atmospheric neutrino masses 
would require a right-handed neutrino with a mass below the GUT scale. 

The see-saw mechanism can be formally derived from the following Lagrangian 

£ = -ULfriLRi'R - ]^VrMrrur + H.c. (3.6) 

where represents left-handed neutrino fields (arising from electroweak doublets), vr rep- 
resents right-handed neutrino fields (arising from electroweak singlets), in a matrix notation 
where the miR matrix elements are typically of order the charged lepton masses, while the 
Mrr matrix elements may be much larger than the electroweak scale, and maybe up to 
the Planck scale. The number of right-handed neutrinos is not fixed, but the number of 
left-handed neutrinos is equal to three. Below the mass scale of the right-handed neutrinos 
we can integrate them out using the equations of motion 

^ = (3.7) 



which gives 

Substituting back into the original Lagrangian we find 



^ - -VLrriLRM^l,, UR = -M^^itiIruI (3.8) 



/: = -hLmLLi^l + H.c (3.9) 



with niLL as in Eq.l 
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3.2 Type II See-Saw and Double See-Saw 



The version of the see-saw mechanism discussed so far is sometimes called the Type I see- 
saw mechanism. It is the simplest version of the see-saw mechanism, and can be thought 
of as resulting from integrating out heavy right-handed neutrinos to produce the effective 
dimension 5 neutrino mass operator 

—HuL^kHuL (3.10) 

where 

K = Yl^M^l^Yi:/ . (3.11) 
One might wonder if it is possible to simply write down an operator by hand similar to 



Eqj3J_0| without worrying about its origin. In fact, historically, such an operator was 
introduced suppressed by the Planck scale (rather than the right-handed neutrino mass 
scales) by Weinberg in order to account for small neutrino masses |3^. The problem 
is that such a Planck scale suppressed operator would lead to neutrino masses of order 
10~^eV which are too small to account for m2 or (though they could account for mi). 
To account for requires dimension 5 operators suppressed by a mass scale of order 
3 X 10^^ GeV if the dimensionless coupling of the operator is of order unity, and the Higgs 
vev is equal to that of the Standard Model. 

One might also wonder if the see-saw mechanism with right-handed neutrinos is the 



only possibility? In fact it is possible to generate the dimension 5 operator in Eq j3.10| by the 
exchange of heavy Higgs triplets of SU{2)l, referred to as the type II see=saw mechanism. 

Alternatively the see-saw can be implemented in a two-stage process by introducing 
additional neutrino singlets beyond the three right-handed neutrinos that we have consid- 
ered so far. It is useful to distingush between "right-handed neutrinos" vr which carry 
B — L and perhaps form SU{2)ii doublets with right-handed charged leptons, and "neutrino 
singlets" S which have no Yukawa couplings to the left-handed neutrinos, but which may 
couple to vr. If the singlets have Majorana masses M55, but the right-handed neutrinos 
have a zero Majorana mass Mrr = 0, the see-saw mechanism may proceed via mass cou- 
plings of singlets to right-handed neutrinos Msr. In the basis (i^^, vr, S) the mass matrix 
is 

/ mrn \ 



rriLR 
niLR Mrs 
\ Mis Mss) 



(3.12) 



Assuming Mss ^ Mrs the light physical left-handed Majorana neutrino masses are then 
doubly suppressed, 

ruLL = niLRM-lMssMls'^ml^ (3.13) 



This is called the double see-saw mechanism. It is often used in string inspired neutrino 
mass models [39|. 
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4. Right-Handed Neutrino Dominance 



In this section we discuss an elegant and natural way of accounting for a neutrino mass 
hierarchy and two large mixing angles, by simply assuming that not all of the right-handed 
neutrinos contribute equally to physical neutrino masses in the see-saw mechanism. This 
mechanism, called sequential dominance, is a technique rather than a model, and can be 
applied to large classes of models. Indeed the conditions for sequential dominance can only 
be understood within particular models, and provide useful clues to the nature of such 
models. 

4.1 Single Right-Handed Neutrino Dominance 

With three left-handed neutrinos and three right-handed neutrinos the Dirac masses tjilr 
are a 3 x 3 (complex) matrix and the heavy Majorana masses M/j/j form a separate 3x3 
(complex symmetric) matrix. The light effective Majorana masses ttill are also a 3 x 3 



(complex symmetric) matrix and continue to be given from Eq. 3^ which is now interpreted 
as a matrix product. From a model building perspective the fundamental parameters 
which must be input into the see-saw mechanism are the Dirac mass matrix niLn and the 
heavy right-handed neutrino Majorana mass matrix Mpm. The light effective left-handed 
Majorana mass matrix miL arises as an output according to the see-saw formula in Eq. |3.5| . 
The goal of see-saw model building is therefore to choose input see-saw matrices mm and 
Mrr that will give rise to one of the successful matrices ttill in Table |l|. 

We now show how the input see-saw matrices can be simply chosen to give the Type 
lA matrix, with the property of a naturally small sub-determinant in Eq. 2.33| using a 



mechanism first suggested in [^]. ^ The idea was developed in |42| where it was called 
single right-handed neutrino dominance (SRHND) . SRHND was first successfully applied 



to the LMA MSW solution in [|3|. 



To understand the basic idea of dominance, it is instructive to begin by discussing a 
simple 2x2 example, where we have in mind applying this to the atmospheric mixing in 
the 23 sector: 

— (;!) <") 

The see-saw formula in Eq j3.5| ttill = ^lrMj^^tuJ^j^ gives: 



=2 



ef , be 



y "I" X y "I" X 











Hi 









mLL= L^i ^ )«( J, h 1 (4.2) 



where the approximation in Eq. 4.2 assumes that the right-handed neutrino of mass Y is 



sufficiently light that it dominates in the see-saw mechanism: 



Y X ' ^ ■ ^ 



The neutrino mass spectrum from Eq.4.2 then consists of one neutrino with mass ~ 
(e^ + f'^)/Y and one naturally light neutrino m2 ^ m^, since the determinant of Eq. |4.2| is 

^See also Q 
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clearly approximately vanishing, due to the dominance assumption [40|. The atmospheric 
angle from Eq.4.2 is tan023 ~ e// which can be large or maximal providing e ~ /, 



even in the case e,f,b <C c that the neutrino Dirac mixing angles arising from Eq.4.1 
are small. Thus two crucial features, namely a neutrino mass hierarchy ^ 7712 and a 
large neutrino mixing angle tan 023 ^ 1; can arise naturally from the see-saw mechanism 
assuming the dominance of a single right-handed neutrino. It was also realised that small 
perturbations from the sub-dominant right-handed neutrinos can then lead to a small solar 



neutrino mass splitting |40|, as we now discuss. 



4.2 Sequential Right-Handed Neutrino Dominance 

In order to account for the solar and other mixing angles, we must generalise the above 
discussion to the 3x3 case. The SRHND mechanism is most simply described assuming 
three right-handed neutrinos in the basis where the right-handed neutrino mass matrix is 



diagonal although it can also be developed in other bases |42, 43]. In this basis we write 
the input see-saw matrices as 

y \ 



Mrr 



(4.4) 



rriLR 



(4.5) 



0X0 
X' J 

^ d a a' 
ebb' 

In [^] it was suggested that one of the right-handed neutrinos may dominante the con- 
tribution to rriLL if it is lighter than the other right-handed neutrinos. The dominance 
condition was subsequently generalised to include other cases where the right-handed neu- 
trino may be heavier than the other right-handed neutrinos but dominates due to its larger 



Dirac mass couplings [42|. In any case the dominant right-handed neutrino may be taken 
to be the one with mass Y without loss of generality. 

It was subsequently shown how to account for the LMA MSW solution with a large 



solar angle |43| by careful consideration of the sub-dominant contributions. One of the 
examples considered in [E^] is when the right-handed neutrinos dominate sequentially, 



\e%\f%\ef\ 
Y 



\xy\ 
X 



> 



\x'y'\ 
X' 



(4.6) 



which is the straightforward generalisation of Eq.4.3 where x,y £ a,b,c and x', y' G a', b' , d . 
Assuming SRHND with sequential sub-dominance as in Eq. 4.f:, then Eqs.3.5, 4.4, give 



g?_ j_ d?_ abide ac , df X 

x"|"y x"|"y X ^ Y » 

6^ I be I ef 

X ^ Y X ^ Y 

• X+^J 



(4.7) 



where the contribution from the right-handed neutrino of mass X' may be neglected ac- 



cording to Eq,4.6. If the couplings satisfy the sequential dominance condition in Eq.4.6 
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then the matrix in Eq.O resembles the Type lA matrix, and furthermore has a naturally 
small sub-determinant as in Eq.|3|. This leads to a full neutrino mass hierarchy 



(4.8) 



and, ignoring phases, the solar angle only depends on the sub-dominant couplings and is 
given by tan ^12 ~ a/(c236 — S23C) |^3|. The simple requirement for large solar angle is then 
a ~ 6- c 

Including phases the neutrino masses are given to leading order in 1712/ tti^ by diago- 
nalising the mass matrix in Eq.4/7 using the analytic proceedure described in Appendix D 
|26]. In the case that d = 0, corresponding to a 11 texture zero in Eq.4.5, we have 

mi ~ O(^) (4.9) 

ma w (4.10) 



XSI2 



+ 1/1^ 



Y 



(4.11) 



where S12 = sin 612 is given below. Note that with SD each neutrino mass is generated by a 
separate right-handed neutrino, and the sequential dominance condition naturally results 
in a neutrino mass hierarchy rrii <^ m2 <^ m-^. The neutrino mixing angles are given to 
leading order in m2/m-i by, 



tan 6*23 
tan 6*12 



\e\ 

m 



^3 



C23\b\ cos{(ph) - S23\c\ cos((;^c) 
,.(<^+</.a-<Ae) Ia|(e*^ + /*c) Y 



(4.12) 
(4.13) 
(4.14) 



[|e|2 + |/|2]3/2X 

where we have written some (but not all) complex Yukawa couplings clS X — I X |e*<P-. The 
phase S is fixed to give a real angle 612 by, 

C23|6| sm{4>b) « S23|c| sin(0c) (4.15) 

where 

4>b = 4>b - 4>a - 4> + S, 

4>c = ^c-^a + (l)e-(l)f-4> + S (4.16) 

The phase (p is fixed to give a real angle ^13 by, 

(f> ^ (pe - (t>a - •/'COSMO (4-17) 

where 

0COSMO = arg(e*6 + /* c). (4.18) 

is the leptogenesis phase corresponding to the interference diagram involving the light- 
est and next-to- lightest right-handed neutrinos pi]. 



- 22 - 



4.3 Types of Sequential Dominance 



Assuming sequential dominance, there is still an ambiguity regarding the mass ordering of 
the heavy Majorana right-handed neutrinos. So far we have assumed that the dominant 
right-handed neutrino of mass Y is dominant because it is the lightest one. We emphasise 
that this need not be the case. The neutrino of mass Y could be dominant even if it is 
the heaviest right-handed neutrino, providing its Yukawa couplings are strong enough to 



overcome its heaviness and satisfy the condition in Eq.4.6. In hierarchical mass matrix 
models, it is natural to order the right-handed neutrinos so that the heaviest right-handed 
neutrino is the third one, the intermediate right-handed neutrino is the second one, and 
the lightest right-handed neutrino is the first one. It is also natural to assume that the 33 
Yukawa coupling is of order unity, due to the large top quark mass. It is therefore possible 
that the dominant right-handed neutrino is the heaviest (called heavy sequential dominance 
or HSD), the lightest (called light sequential dominance or LSD), or the intermediate one 
(called intermediate sequential dominance or ISD). This leads to the six possible types of 
sequential dominance corresponding to the six possible mass orderings of the right-handed 
neutrinos as shown in Table|l|. In each case the dominant right-handed neutrino is the one 
with mass Y, and the leading subdominant right-handed neutrino is the one with mass 
X. The resulting see-saw matrix miL is invariant under re-orderings of the right-handed 
neutrino columns, but the leading order form of the neutrino Yukawa matrix 1^ is not. 

It is worth emphasising that since all the forms above give the same light effective see- 
saw neutrino matrix rnxL in Eq j4.7l , under the sequential dominance assumption in Eq.4.6 



this implies that the analytic results for neutrino masses and mixing angles applies to all 
of these forms. They are distinguished theoretically by different preferred leading order 
forms of the neutrino Yukawa matrix shown in the table. These leading order forms 
follow from the the large mixing angle requirements e ~ / and a ~ 6 — c. ^ Thus we see 
that LSDa, and ISDa are consistent with a form of Yukawa matrix with small Dirac mixing 
angles, while HSDa and HSDb correspond to the so called "lop-sided" forms. LSDb and 
ISDb correspond to the D-brane inspired "single right-handed democracy" form studied in 



1 44]. They are also distinguished by leptogenesis and lepton flavour violation as we shall 
see. 

For example, suppose that we impose the theoretical requirement that the neutrino 
Yukawa matrix resembles hierarchical quark matrices, and have a large 33 element of order 
unity, but no other large off-diagonal entries. Then the large mixing angle requirements 
e ~ / and a ~ 6 — c immediately excludes HSDa, HSDb, LSDb and ISDb. We are left 
with LSDa and ISDa as the remaining possibilities. If we further impose the requirement 



of a 11 texture zero, as motivated by the GST relation |45|, then a ~ 6 — c excludes ISDa, 
and we are left uniquely with LSDa. We shall later discuss an example of a realistic model 
of all quark and lepton masses and mixing angles based on LSDa. For now we note that 



that for LSDa in order to satisfy the sequential dominance condition in Eq.4.6 the heavy 



^Note that the leading order in the Table only gives the independent order unity entries in the matrix, 
so that for example in LSDb we would expect 6 — c ~ 1 in general, and not zero. 
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Table 2: Types of sequential dominance (SD), classified according to the mass ordering of the right- 
handed neutrinos. Light sequential dominance (LSD) corresponds to the dominant right-handed 
neutrino of mass Y being the lightest. Intermediate sequential dominance (ISD) corresponds to the 
dominant right-handed neutrino of mass Y being the intermediate one. Heavy sequential dominance 
(HSD) corresponds to the dominant right-handed neutrino of mass Y being the heaviest. 
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Majorana masses must be necessarily strongly hierarchical, 

y < X < X'. (4.19) 

The reason is that the heavy right-handed neutrino of mass X' has order unity Yukawa 
couplings to left-handed neutrinos, which implies that the lightest right-handed neutrino 
of mass Y must be significantly lighter in order to dominate. 

4.4 Leptogenesis Link 

It is interesting to note that in LSDa, assuming a 11 texture zero, there is a link between 
the CP violation required for leptogenesis, and the phase 6 measurable in accurate neutrino 
oscillation experiments. This can be seen from Eq. 4.18| which may be expressed as 



+ J. ~ \b\s23S2 + c C23S3 

tan (/)coSMO ~ 77] m • (4-20) 

\0\S23C2 + C C23C3 
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From Eqs, 4.17, 4.15, 4. 16 



tan((/)cosMo + <5) ~ —77] nn (4-21) 

-|0|C23C2 + |c|S23C3 

where we have written Si = sin r]i,Ci = cos rji where 

7/2 = 4- 4>e, V3 = 4>c- 4>f (4.22) 

are invariant under a charged lepton phase transformation. The reason that the see-saw 
parameters only involve two invariant phases r]2 , % rather than the usual six is due to the 
LSD assumption which has the effect of decoupling the heaviest right-handed neutrino, 
which removes three phases, together with the assumption of a 11 texture zero, which 
removes another phase. 



Eq. 4.21 shows that (5 is a function of the two see-saw phases r/2 , % that also determine 



'/'COSMO ill Eq. 4.20| . If both the phases are zero, then both (/^COSMO and 5 are 



necessarily zero. This feature is absolutely crucial. It means that, barring cancellations, 
measurement of a non-zero value for the phase J at a neutrino factory will be a signal of a 
non-zero value of the leptogenesis phase ^COSMO- We also find the remarkable result 

I '/'COSMO I = \(Pl3l30u\- (4.23) 

where 4'i3f30u is the phase which enters the rate for neutrinoless double beta decay |p6|| . 

4.5 Comparison to the Smirnov Approach 

An early approach to obtaining large mixing angles from the see-saw mechanism was pro- 
posed by 1 47 1, which is sometimes confused with sequential dominance. The purpose of 



this subsection is to briefly review the Smirnov approach, and explain how it differs from 
sequential dominance. The Smirnov approach for obtaining large mixing angles from the 
see-saw mechanism, is based on the theoretical assumption of having no large mixing angles 



in the Yukawa sector [47| 



We shall briefly discuss the two family case considered in |47]. For this case the physical 
lepton mixing angle is written as 

e = ei^^ - ei^^ + Oss (4.24) 

where ^2^^ is the left-handed mixing angle which diagonalises the charge lepton Yukawa 
matrix, O^j^^ is the left-handed mixing angle which diagonalises the neutrino Yukawa ma- 
trix, and 6ss is defined to be the additional angle which results from the presence of the 



see-saw mechanism. The basic idea [47| was that a large mixing angle 6 could originate 
from the see-saw mechanism via Oss with ^2^^ and ^^^^ being small. ^ Smirnov obtains 
an approximate analytic expression for Oss in the two family case. 



tan(0^^« - 0^««) 
tan2(0^^« -0^««) +eA^ 



tan 055 ^ -26^ .„„2..^.." a^nn^ 'm (4-25) 



*Note that this is not a requirement of sequential dominance, ahhough it may be satisfied by LSDa or 
ISDa, as discussed previously. 
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where O'j^^ is the mixing angle which diagonahses the neutrino Yukawa Y^i? matrix on the 
right, is the mixing angle which diagonahses the heavy Majorana neutrino matrix 

Mrr, is the ratio of neutrino Yukawa (Dirac) matrix eigenvalues, and is the ratio 
of heavy Majorana matrix eigenvalues. The conditions that 6ss is large are 

tan(^^^« < (4.26) 
e^^ < {e^f (4.27) 

which, for a typical quark- like hierarchy ^ 1, implies both a very accurate equality of 
mixing angles O'jf'^ = O^r^ and very strongly hierarchical heavy Majorana masses (much 
stronger than the Dirac mass hierarchy). 

The conditions in Eqs. 4.26| , |4.27] are clearly nothing to do with sequential dominance in 



general. For one thing since some versions of sequential dominance involve large neutrino 
Yukawa mixing angles ^^^^ and do not require Oss to be large, which is the basic assumption 
of this approach. However there are classes of sequential dominance model such as LSDa 
where O^j^^ is small and Ogs is large. Furthermore in this class of model there is a strong 
hierarchy of Majorana masses. One might be tempted to think that LSDa is the same 
as the Smirnov approach, and this has led to some confusion in the literature which we 



would like to clear up here. The important point to emphasise is that |47] never talks 
about one of the right-handed neutrinos giving the dominant contribution to the heaviest 
physical neutrino via the see-saw mechanism, or indeed about the relative contribution of 
the right-handed neutrinos to the see-saw mechanism in general. Thus there is no natural 
mechanism present for generating a neutrino mass hierarchy in [47|, which is concerned 



only with the condition for generating large mixing angles. The point about sequential 
dominance is that it can naturally generate a neutrino mass hierarchy and large mixing 
angles, as simple ratios of Yukawa couplings of dominant and subdominant right-handed 
neutrinos. 

A simple counter example will illustrate this point. Condider the following matrices, 

where Aij,aij are order unity coefficients. These matrices clearly satisfy the conditions 
in Eqs j4.26 , 4.27 , since 9'^^ ~ 9'^^ ^ €d and €m ~ However these matrices do not 



satisfy the dominance conditions. Both right-handed neutrinos will contribute equally at 
0(1/M) via the see-saw mechanism to the heaviest physical neutrino mass. Without the 
dominance of a single right-handed neutrino the neutrino mass hierarchy will require some 
tuning. The tuning required for the atmospheric mixing angle involving second and third 
families with be rather mild since 7712 /ma is not so small, however when this scheme is 
extended to all three families, further tuning will be required to obtain a large solar mixing 
angle in a natural way. In actual examples given in even more tuning is likely to be 
required since the angles 9'^^, 9'^^ were both independently supposed to be larger than 
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The conclusion is that Smirnov's approach did not recognise right-handed neutrino 
dominance, contrary to some recent claims in the literature, but it does provide a com- 
plementary approach to a large mixing angles from the see-saw mechanism. At first sight 
it appears to have some similarities to LSDa, however without the missing ingredient of 
sequential dominance, to achieve two large mixing angles together with a neutrino mass 
hierarchy will require some degree of fine-tuning. The conditions proposed by Smirnov are 
therefore neither necessary nor sufficient for right-handed neutrino dominance. 



5. See-Saw Standard Models 

In this section we show how the see-saw mechanism can be accomodated in the Standard 
Model and its Supersymmetric Extension, where it leads to lepton flavour violation. 

5.1 Minimal See-Saw Standard Model 

We now briefly discuss what the standard model looks like, assuming a minimal see-saw 
extension. In the standard model Dirac mass terms for charged leptons and quarks are 
generated from Yukawa couplings to Higgs doublets whose vacuum expectation value gives 
the Dirac mass term. Neutrino masses are zero in the Standard Model because right- 



handed neutrinos are not present, and also because the Major ana mass terms in Eq.|3.1 
require Higgs triplets in order to be generated at the renormalisable level. The simplest 
way to generate neutrino masses from a renormalisable theory is to introduce right-handed 
neutrinos, as in the Type I see-saw mechanism, which we assume here. The Lagrangian 
for the lepton sector of the standard model containing three right-handed neutrinos with 
heavy Major ana masses is ^ 



+ H.C. (5.1) 



where eab = — Cba, ei2 = l, and the remaining notation is standard except that the 3 right- 
handed neutrinos have been replaced by their CP conjugates f ?, and is a complex 
symmetric Majorana matrix. When the two Higgs doublets get their VEVS < i^^ >= ^2-, 
< H\ >= f 1, where the ratio of VEVs is defined to be tan /3 = ^2/^1, we find the terms 

Cmass = -vi^eie] - V2^v,v''^ - \MtRvtv''^ + H.c. (5.2) 

Replacing CP conjugate fields we can write in a matrix notation 

C-mass = -eLVlY'^*eR - l'LV2Y''*l'R - ]^v]^M*^rVr + H.c. (5.3) 

It is convenient to work in the diagonal charged lepton basis 

diag(me,m^,nv) = V^^i'i^'^VJ^ (5.4) 

^We introduce two higgs doublets to pave the way for the supersymmetric standard model. For the same 
reason we express the standard model Lagrangian in terms of left-handed fields, replacing right-handed fields 
by their CP conjugates. In the case of the minimal standard see-saw model with one Higgs doublet one of 
the two Higgs doublets by the charge conjugate of the other, Hd = H'^. 
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and the diagonal right-handed neutrino basis 

diag(Mi, M2, M3) = V,„M*rrVJ^ (5.5) 

where VeL,VeR,Vyj^ are unitary transformations. In this basis the neutrino Yukawa cou- 
phngs are given by 

Y- = Ve^Y^'vJ^ (5.6) 

and the Lagrangian in this basis is 

C-mass =- {eLp.LfL)diag{me,m^,mr){eRfJ.RTR)'^ 

- {^eLVf,LVTL)yV2{vRlVR2VR2f 

- {iyRliyR2l^R3)diag{Mi,M2,M3){uRiUR2UR3)'^ + H.C. (5.7) 

After integrating out the right-handed neutrinos (the see-saw mechanism) we find 
Cmass = - (eL//LrL)diag(me,m^,nv)(eR//i?Tij)^ 

- ^{l'eLl't^Ll^TL)mLL{T^ell'fill^TlV + H.C. (5.8) 

where the hght effective left-handed Majorana neutrino mass matrix in the above basis is 



given by the following see-saw formula which is equivalent to Eq.3.5, 

rriLL = -vlVdiagiU^^ ,M^^ ,M^^)Y''^ (5.9) 
Eq j5.8| is equivalent to Eq. 2.22 when expressed in the charged lepton mass basis, which we 



have derived starting from the standard model Lagrangian using the see-saw mechanism. 
5.2 Minimal Supersymmetric See-Saw Standard Model 

It is well known that large mass scales such as are required in the see-saw mechanism can 
be stabilised by assuming a TeV scale N=l supersymmetry which cancels the quadratic 
divergences of the Higgs mass. Thus it is natural to generalise the see-saw standard model 
to include supersymmety. When this is done the leptonic part of the superpotential with 
three right-handed neutrinos is given by 

where €ab = —^ba and ei2 = 1. The SU{2) representations of the lepton superfield doublets 
can be expressed as follows (suppressing family indices for simplicity): 



U=\.'\. (5.11) 




The superfields are defined in the standard way as follows (suppressing gauge indices): 

= (z>l,,i/lJ 
Ci = (eL,,eLj 
ef = (ei,,eij 



Hd = {Hd,Hd), (5.12) 
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with i,j = 1 ... 3 labeling family indices. The soft breaking Lagrangian Cgoft in the lepton 
sector takes the form (dropping "helicity" indices): 



1 



(5.13) 



soft ^aoi^^d^i^ '^i]^] I ^■^u^i-'H]'^ ] 

The Yukawa terms in the Lagrangian are given from the superpotential by replacing 
two of the superfields by their fermion components, and one of the superfields by its 
scalar component, and including an overall minus sign. Then the leptonic part of the 



superpotential in Eq.5.1C reduces to the standard model lagrangian in Eq.5.1, and the 
discussion then follows that of the previous section. For the charged leptons, we have as 
before 



in which 



diag (me, m/„m^) = VelvIY*VI^] 
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(5.14) 



(5.15) 



The important new feature provided by SUSY is the existence of scalar partners to 
the leptons (sleptons) which can give lepton flavour violating (LFV) effects, which arise 
as discussed in the following. To discuss these effects we first need to express the sleptons 
in terms of their mass eigenstates. It is usually convenient however to begin by rotating 
the sleptons in exactly the same way as the lepton. In this basis, which we call the MNS 
basis, the photino interactions conserve flavor, while the wino (and higgsino) interactions 
violate flavor by U, in analogy to the gauge boson interactions in the SM. Therefore, the 
diagonalization of the scalar mass matrices proceeds in two steps. First, the sleptons are 
rotated "parallel" to their fermionic superpartners; i.e., we do unto sleptons as we do unto 
leptons: 
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= Ve^ 
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(5.16) 



where in the MNS basis the slepton fields {cl, p-L,?!,) are SUSY partners of the physical 
mass eigenstate quarks (eL,^L,rL), respectively, (i.e. (ei,,eL) share the same superfield 
where both components of the superfield have been subject to the same rotation, thereby 
preserving the superfield structure), and similarly for the other terms. 

The slepton fields expressed in the MNS basis are often more convenient to work with, 
even though they are not mass eigenstates. Their 6x6 mass matrices are obtained by 
adding the electroweak symmetry breaking contributions and then rotating to the MNS 
basis. They have the following form: 



nif 



^t>LR 



cos 2/3 



(M| + 2M2 )i 



K)l 



tan (3iJ,*me 



("i|)lr 
iml)RR + ml- 



tan PfiTJie 



cos 2/3 



Mi sin^ 
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in which 6w is the electroweak mixing angle, 1 stands for the 3x3 unit matrix, and we 
have written me = diag^mg, m^,mr)- The flavor-changing entries responsible for lepton 
flavour violation are contained in the off-diagonal entries of the soft slepton mass matrices 
above, which are given by 

iml)LL = VE.mfvl^ {ml)RR = VE^m%V^^ = vIVeJ^'V^^. (5.18) 

5.3 Lepton Flavour Violation 

The renormalisation group equations (RGEs) contain additional terms relative to the 
MSSM. The additional terms imply that even if the soft slepton masses are diagonal at the 
GUT scale, then in this case we would find that three separate lepton numbers Le,L^,Lr 
are not conserved at low energies, since the new RGE terms do not preserve these symme- 
tries in general if there are right-handed neutrinos below the GUT scale. Below the mass 
scale of the right-handed neutrinos we must decouple the heavy right-handed neutrinos 
from the RGEs, and then the RGEs return to those of the MSSM. Thus the lepton number 
violating additional terms are only effective in the region between the GUT scale and the 
mass scale of the lightest right-handed neutrino, and all the effects of lepton number vio- 
lation are generated by RGE effects over this range. The effect of RGE running over this 
range will lead to off-diagonal slepton masses at high energy, which result in off-diagonal 
slepton masses at low energy, and hence observable lepton flavor violation in experiments. 

Assuming universal soft parameters at Mgut-, "ii(0) = m^(0) = m^I, where / is the 
unit matrix, and ^,^(0) = AY^, the renormalisation group equation (RGE) for the soft 
slepton doublet mass may be written as 



dml _ f dm\ \ (3m2 + A^) 

dt ~ V dt )y=o 16vi^ 



YuY^] (5.19) 



where in the basis in which the charged lepton Yukawa couplings are diagonal, the first term 
on the right-hand side is diagonal. In running the RGEs between Mgut and a right-handed 
neutrino mass Mj the neutrino Yukawa couplings lead to an approximate contribution to 
the slepton mass squared matrix of 

" "1^ (%^) + [^^^A 

This shows that, to leading log approximation, off-diagonal slepton masses may be gen- 
erated depending on the form of the neutrino Yukawa matrix. The off-diagonal slepton 
masses give rise to LEV processes such as n ^ e'j, r — > jij, r — > 67. From a future 
observation of these processes one may infer information about the slepton mass matrix, 
and then use this information to make inferences about the neutrino Yukawa matrix, and 
heavy right-handed neutrino masses. This proceedure would be impossible in the SM, 
and is an example of how supersymmetry may in the future provide a window into the 
Yukawa matrices which would not otherwise be possible. This was originally discussed in 
1 48, 49, 51 1 and has been discussed recently in |5^, ^ 55, 56, 57, 58 1. 
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At leading order in a mass insertion approximation the branching fractions of LFV 
processes are given by 

3 

BR{li ^ /j7) « ^/(M2,/i,ms)|m? |2tan2/5 (5.21) 

Up 

where li = e,l2 = ^3 = t, and where the off-diagonal slepton doublet mass squared is 
given in the leading log approximation (LLA) by 

2(LLA) ^ _ (3mg+^g) 
Lij Svr^ 



where in sequential dominance, in the notation of Eqs,4.4,4.5 the leading log coefficients 



relevant for /i — > 67 and r fij are given approximately as 

C21 = ao In — — + de In —— 
A Y 

C32 = be In — + ef In — (5.23) 

A global analysis of LFV has been performed in the constrained minimal supersym- 
metric standard model (CMSSM) for the case of sequential dominance, focussing on the 
two cases of HSD and LSD [^]. The results for HSD show a large rate for r which 



is the characteristic expectation of lop-sided models in general |53] and HSD in particular. 



The results are based on an exact calculation, and the error incurred if the LLA were used 



can be as much as 100% |57]. For LSD r — > is well below observable values. Therefore 
T ^ provides a good discriminator between the HSD and LSD types of dominance. 
The rate for fj, ^ can be large or small in each case. 

6. GUTs and Family Symmetry 

We have seen that atmospheric neutrino masses would seem to require a right-handed 
neutrino with a mass below the GUT scale. Such a mass scale demands an explanation, 
and in fact one must then explain why the right-handed neutrinos are so light compared 
to the Planck scale. In order to explain this, one clearly needs a theory of right-handed 
neutrino masses capable of protecting the right-handed neutrino masses by some symmetry 
which is subsequently broken at some scale. Suitable symmetries can correspond to either 
unification or family symmetries, as we now discuss. 

6.1 Models Based on GUTs and Family Symmetry 

One of the exciting things about the discovery of neutrino masses and mixing angles is 
that this provides additional information about the flavour problem - the problem of un- 
derstanding the origin of three families of quarks and leptons and their masses and mixing 
angles (Fig.^). Early approaches to the problem of quark masses and mixing angles in- 
cluded the postulate that some entries in the Yukawa matrices were equal to zero (the 



so-called "texture zeroes") thereby reducing the number of free parameters [59|. In this 



approach the quark and lepton Yukawa matrices are assumed to be hierarchical in nature 
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with an order unity entry in the 33 entry. Another complementary approach is to assume 
that the Yukawa matrices are democratic with order unity entry everywhere [Bol , and both 



approaches have been followed for neutrino masses and mixings |6^, 34, A 
specific model of the neutrino mass matrix with texture zeroes, but with a texture zero 
in the 33 position was proposed by Zee [^], and this has been developed recently by a 
number of authors |7^. Unfortunately the simplest Zee texture is now excluded 

by experiment, although a non-minimal Zee type model remains viable |fl| . 

To understand the origin of the postulated forms of Yukawa matrices, one must appeal 
to some sort of Family symmetry Gpamiiyj which acts in the direction shown in Fig.Q. In the 
framework of the see-saw mechanism, new physics beyond the standard model is required 
to violate lepton number and generate right-handed neutrino masses which are typically 
around the GUT scale. This is also exciting since it implies that the origin of neutrino 
masses is also related to some GUT symmetry group GguT) which unifies the fermions 
within each family as shown in Fig.^. 

Putting these ideas together we are suggestively led to a framework of new physics 
beyond the standard model based on N=l SUSY ^ with commuting GUT and Family 
symmetry groups, 

CgUT X GpAM (6.1) 

There are many possible candidate GUT and Family symmetry groups some of which 
are listed in Table |3[ Unfortunately the model dependence does not end there, since the 
details of the symmetry breaking vacuum plays a crucial role in specifying the model and 



determining the masses and mixing angles, resulting in many models as given in [72| - |112| 
(listed alphabetically). These models may be classified according to the particular GUT 
and Family symmetry they assume as shown in Table |3[ 
We have used the notation that 

51 = SU{5) X U{1) (6.2) 

422 = SU{A)ps X SU{2)l x SU{2)r (6.3) 

3221 = SU{3)c X SU{2)l x SU{2)r x U{1)b-l (6.4) 

321 = SU{3)c X SU{2)l x U{1)y (6.5) 

where 422 is the Pati-Salam gauge group, 3221 is the left-right symmetric gauge extension, 
321 is the Standard Model gauge group. 

Another complication is that the masses and mixing angles determined in some high 
energy theory must be run down to low energies using the renormalisation group equations 



(RGEs) |113| , |ll4 |117| , |115| , |119|] . Large radiative corrections are seen when the see-saw 
parameters |104] are tuned, since the spectrum is sensitive to small changes in the pa- 
rameters, and this effect is sometimes used to magnify small mixing angles into large ones 
I lis, 35, 116, |118(| . This idea has however been criticised in [ [12[| . In natural models based 



on SRHND the parameters are not tuned, since the hierarchy and large atmospheric and 



''Supersymmetry enables the gauge couplings to meet at the GUT scale to give a self-consistent unifica- 
tion picture. 
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Table 3: Some candidate GUT and Family symmetry groups, and the papers that use these 
symmetries to successfully describe the LMA MSW solar solution and the atmospheric neutrino 
data. 
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solar angles arise naturally as discussed in the previous section. Therefore in SRHND 
models the radiative corrections to neutrino masses and mixing angles are only expected 



to be a few per cent, and this has been verified numerically [121|. 



6.2 50(10) X U{1) 

As an example we shall here consider a model based on a GUT group 50(10) and a family 
symmetry U{1). We shall suppose that the GUT symmetry is broken via a Pati-Salam 



group and define the model in terms of the subgroup 422 [94|. This model provides an 
example of the use of both the ^7(1) family symmetry to generate inter- family hierarchies, 
and the use of Clebsch-Gordon coefficients from the GUT group to generate intra-family 
structure. 

The left-handed quarks and leptons are accommodated in the following 422 represen- 
tations, 

\ U U U V j 

where i = 1 ... 3 is a family index. The Higgs fields are contained in the following repre- 
sentations, 

(where hi and h2 are the low energy Higgs superfields associated with the MSSM.) 
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The quark and lepton mass hierarchies 



Charge +2/3 quarks 



Charge -1/3 quarks 



Charged leptons 



A neutrino 
hierarchy 




Figure 4: The fermion masses are here represented by a lego plot. We have multiplied the masses 
of the bottom, charm and tau by 10, the strange and muon by 10^, the up and down by 10^, the 
electron by 10^ to make the lego blocks visible. It is natural to assume a normal neutrino hierarchy. 
We have multiplied the third neutrino mass by 10^^ and the second neutrino mass by 10^^ to make 
the lego blocks visible. This underlines how incredibly light the neutrinos are. The symmetry 
groups Ggut and Gpamiiy act in the directions indicated. 



The two heavy Higgs representations are 
H^' = (4,1,2) = I 



'^H '^H '^H 
'^H '^'h 



and 



U§ U% VH 



(4,1,2): 

The Higgs fields are assumed to develop VEVs, 

< H >=< UH >~ Mgut, < H >=< VH >~ Mgut 
leading to the symmetry breaking at Mgut 

SU(4) O SU(2)^ O SU(2)^ — > SU(3)p SU(2)^ ® U(l)y 



(6.9) 



(6.10) 



(6.11) 



(6.12) 



in the usual notation. Under the symmetry breaking in Eq j6.12 , the Higgs field h in Eq 
splits into two Higgs doublets hi, h2 whose neutral components subsequently develop weak 
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scale VEVs, 

<,h1>=vi, <hl>=V2 (6.13) 

with tan /3 = vilv\. 

To construct the quark and lepton mass matrices we make use of non-renormahsable 



operators [122] of the form: 



-w) (m) 

The 9 fields are Pati-Salam singlets which carry U{1) family charge and develop VEVs 
which break the U{i) family symmetry. They are required to be present in the operators 
above to balance the charge of the invariant operators. After the H and 9 fields acquire 
VEVs, they generate a hierarchy in i) effective Yukawa couplings and ii) Majorana masses. 
These operators are assumed to originate from additional interactions at the scale M > 
Mgut- The value of the powers pij and qij are determined by the assignment of U{1) 
charges, with Xg = -1 then pij = {X^^ + X^j + Xh) and qij = {X^i + X^j + Xh)- 

The contribution to the third family Yukawa coupling is assumed to be only from the 
renormalisable operator with n = p = leading to Yukawa unification. The contribution 
of an operator, with a given power n, to the matrices Yf=u,d,i^,e, ^Irr is determined by 
the relevant Clebsch factors coming from the gauge contractions within that operator. A 



list of Clebsch factors for all n = 1 operators can be found in the appendix of . These 
Clebsch factors give zeros for some matrices and not for others, hence a choice of operators 
can be made such that a large 23 entry can be given to and not Yu^d,e- We shall write, 

then we can identify 5 with mass splitting within generations and e with splitting between 
generations. 

The choice of ^7(1) charges are as in [94| and can be summarised as X^i = (1,0,0), 



X^i = (4,2,0), Xh = 0, Xh = and Xjj = 0. This fixes the powers of e in each entry 
of the Yukawa matrix, but does not specify the complete operator. The Yukawa couplings 



are specified by a particular choice of operators, | 123 , 94] with the property 

O ~ {HH) ~ (5, O' ~ {HHf ~ (5^, O" ~ {HHf ~ 5^. (6.17) 
The Clebsch factors play an important part in determing the form of the Yukawa matrices. 



The particular operator choice in |94] leads to the quark and lepton mass matrices below. 
For example the Clebsch coeficients from the leading operator in the 22 positon gives the 
ratio : 1 : 3 in the Yu^£,^e matrices. This ratio along with subleading corrections provides 



the correct rric ■ rag : ratio ]124] 
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The final form of tlie Yukawa matrices is 112 



ye 









rt 




5^5 



3 
52 £2 25 



(6.18) 



(6.19) 



(6.20) 



(6.21) 





where the numerical Clebsch factors are displayed explicitly, and rt,ri,,ri, are order unity 



parameters which quantify the deviations from exact Yukawa unification 123 |, but all other 
order unity coefficients have been dropped. 

The Majorana operators are assumed to arise from an m = operator in the 33 
position and m = 1 operators elsewhere, resulting in 



Mrr 



f de^ 6e^ 6e^ ' 
6e^ 6e^ 6e^ 
\Se^ 6e^ 1 



M3. 



(6.22) 



In the neutrino sector the matrices above satisfy the condition of sequential dominance 
in which a neutrino mass hierarchy naturally results with the heaviest (third) right-handed 
neutrino being mainly responsible for the atmospheric neutrino mass, and the second heav- 
iest right-handed neutrino being mainly responsible for the solar neutrino mass. Thus this 
model corresponds to HSDa in Table |2|. Using the HSDa ordering in Table || with the 
matrices in Eqs. 3.21 , |6.22 we can use the analytic results in Eqs. 4.9- 4.14 to give estimates 
of neutrino masses 



mi 



m2 



and neutrino mixing angles: 



tan 6*23 



Ms 
i4S^ + rl 



26 



M3 



(6.23) 
(6.24) 
(6.25) 

(6.26) 
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tan 6v 



12 



713 



12 r2^ I 25 



((25)2+r2)3/2 



(6.27) 



(6.28) 



which are a good fit to the LMA MSW solution for e and 6 as in Eg. 3.16 . 
6.3 50(10) X SU{3) 

As an example of a model based on a non-Abelian family symmetry, we briefly review the 
model proposed in |96]. The model uses the largest family symmetry SU{3) consistent 
with 5*0(10) GUTs. An important further motivation for SU{3) family symmetry is, 
in the framework of sequential dominance, to relate the second and third entries of the 
Yukawa matrix, as required to obtain an almost maximal 23 mixing in the atmospheric 
neutrino sector [^]. In this framework we already saw that the theoretical requirements 
that the neutrino Yukawa matrix resembles the quark Yukawa matrices, and therefore has 
a large 33 element with no large off-diagonal elements and a texture zero in the 11 position 
||45[| , leads uniquely to LSDa in Table ^, where the dominant right-handed neutrino is the 
first (lightest) one. Assuming this then the atmospheric neutrino mixing angle is given 
by tan023 ^ ^2^1 ^ 1- "^^^ sequential dominance conditions which were assumed in 
Eq j4.6| will here be derived from the symmetries of the model. Thus this model provides 
an example of the application of sequential dominance to realistic models of flavour, and 
shows how the conditions of sequential dominance which were simply assumed earlier can 
motivate models based on GUTs and family symmetry which are capable of explaining 
these conditions. In other words, the conditions for sequential dominance can provide 
clues to the nature of the underlying flavour theory. 

The starting point of the model is the observation that an excellent fit to all quark 
data is given by the approximately symmetric form of quark Yukawa matrices []65| 

(6.29) 

1 J 

where the expansion parameters e and e are given by 

ewO.05, €^0.15. (6.30) 
This motivates a particular model in which the three families are unified as triplets 




under an SU{3) family symmetry, and 16's under an 50(10) GUT |9^, |107| , |96| , 

V'i = (3,16), (6.31) 
where as before the 50(10) is broken via the Pati-Salam group giving the equivalent 422 

^, = (3,4,2,1), V5i = (3,4,1,2). (6.32) 



reps in Eqs.5.q 6.7], 
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Further symmetries Rx Z2 x f/(l) are assumed to ensure that the vacuum aUgnment leads 
to a universal form of Dirac mass matrices for the neutrinos, charged leptons and quarks 
1 96]. To build a viable model we also need spontaneous breaking of the family symmetry 



SU{3) — > SU{2) — > Nothing 



(6.33) 



To achieve this symmetry breaking additional Higgs fields ^3, 03 , (/>23 and 023 ^i'^ required. 
The largeness of the third family fermion masses implies that SU{3) must be strongly 
broken by new Higgs antitriplet fields 03 which develop a vev in the third SU{3) component 
< 03 >^= (0,0,03) as in 03 transforms under SU{2)ji as 3 © 1 rather than being 

5C/(2)r singlets as assumed in [^], and develops vevs in the SU{3) x SU{2)ji directions 



< 03 > = < 03 > = 




(6.34) 



The symmetry breaking also involves the SU (3) antitriplets 023 which develop vevs |95[| 



< 023 > = 



(6.35) 



where, as in |95], vacuum alignment ensures that the vevs are aligned in the 23 direction. 
Due to D-flatness there must also be accompanying Higgs triplets such as 023 which develop 
vevs Q 



< 023 > = 



/o\ 

1 

v) 



b. 



(6.36) 



We also introduce an adjoint S field which develops vevs in the SU (4) ps x SU {2)ji direction 
which preserves the hypercharge generator Y = T^ji + (B — L)/2, and implies that any 
coupling of the S to a fermion and a messenger such as Ti'^'j^ip^^x^^ ^ where the SU{2)ji and 
SU{A)ps indices have been displayed explicitly, is proportional to the hypercharge Y of 
the particular fermion component of ij)'^ times the vev a. In addition a 9 field is required 
for the construction of Majorana neutrino masses. 

The leading operators allowed by the symmetries are 



Yuk 



1 

M2 
S 

M3 

1 



il;i(j)lijjj(j)lh 

V'j</'23V'j023^ 



(6.37) 
(6.38) 
(6.39) 



where the operator mass scales, generically denoted by M may differ and we have sup- 
pressed couplings of 0(1). 
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The final form of tlie Yukawa matrices and heavy Major ana matrix after inserting a 
particular choice of order unity coefficients is [p6i 



/o 



yd 



ye 



/O 



/o 



/o 



-0.9e^ 

-0.7e3 

-1.6e3 
-0.7e3 

-1.2e3 
-0.9e3 



3^ 
3'' 



-2 



+ e2 



1.6e^ 0.7e'^ 



3e^ 
3e-2 

1.2e' 



3e^ 



-ae 



0.9e3 

— ae^ 
1 



e. 




(6.40) 



(6.41) 



(6.42) 



(6.43) 



(6.44) 



The model gives excellent agreement with the quark and lepton masses and mixing 



angles. For the up and down quarks the form of Y"^ and Y'^ given in Eq. |6.40 , 6.41 is 



consistent with the phenomenological fit in Eq.6.2£. The charged lepton mass matrix is 



of the Georgi-Jarslkog |124] form which, after including radiative corrections, gives an 
excellent description of the charged lepton masses. In the neutrino sector the parameters 
satisfy the conditions of sequential dominance |4.6| , with the lightest right-handed neutrino 
giving the dominant contribution to the heaviest physical neutrino mass, and the second 
right-handed neutrino giving the leading subdominant contribution, providing that a ~ e. 
It thus falls into the category of LSDa in Table |2[ 

Analytic estimates of neutrino masses and mixing angles for sequential dominance were 
derived in [p^], and for the special case here of LSDa, with the 11 neutrino Yukawa coupling 



the neutrino masses are obtained, 



and neutrino mixing angles: 



mi 


~ e-2 


^2 
Ms 




« 5. 


Ms 


1713 


^ 15 


^2 

Ms 


tan 


^23 ^ 


1.3 


tan 


^12 ^ 


0.66 




^13 ^ 


1.6e 



(6.45) 
(6.46) 
(6.47) 

(6.48) 
(6.49) 
(6.50) 
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Note that the physical lepton mixing angle ^13 receives a large contribution from the 
'13 ~ high energy scale, for this choice of parameters, compared 



neutrino sector 



to the current CHOOZ limit ^13 < 0.2 [14|. However the physical mixing angles will receive 
charged lepton contributions |26| and all the parameters are subject to radiative corrections 
in running from the high energy scale to low energies, although in sequential dominance 
models these corrections are only a few per cent |121|. Thus the model predicts that ^13 is 



close to the current CHOOZ limit, and could be observed by the next generation of long 
baseline experiments such as MINOS or OPERA. 



7. Conclusions 



This focus of the review has been on "mainstream" neutrino mass models, defined as see- 
saw models involving three active neutrinos which are capable of describing both the atmo- 
spheric neutrino oscillation data, and the LMA MSW solar solution, which is now uniquely 
specified by recent data. We have briefly reviewed the current experimental status, showed 
how to parametrise and construct the neutrino mixing matrix, and presented the leading 
order neutrino Majorana mass matrices. We then introduced the see-saw mechanism, and 
discussed a natural application of it to current data using the sequential dominance mech- 
anism, which we compared to an early proposal for obtaining large mixing angles. We 
showed how both the Standard Model and the Minimal Supersymmetric Standard Model 
may be extended to incorporate the see-saw mechanism, and showed how the latter case 
leads to the expectation of lepton flavour violation. The see-saw mechanism motivates 
models with additional symmetries such as unification and family symmetry models, and 
we tabulated some possible models, before focussing on two particular models based on 
50(10) grand unification and either U{1) or SU{3) family symmetry as specific examples. 
We have provided extensive appendices which include techniques for analytically diago- 
nalising different types of mass matrices involving two large mixing angles and one small 
mixing angle, to leading order in the small mixing angle. 

Neutrino physics has witnessed a renaissance period with the watershed provided by 
Super-Kamiokande in 1998. Before then we did not know whether atmospheric and solar 
neutrino oscillations were fact or fancy. Now we know they are fact whose explanation 
requires two large mixing angles. We have seen in this review that there is no shortage of 
theoretical models which can account for these data. Even discounting those theoretical 
models which predicted a small solar angle, or vacuum oscillation mass splittings, there 
are many many models that can describe the current data. In this review we have tried 
to emphasise useful approaches and techniques, rather than giving a detailed catalogue of 
all possible models. We make no apology for emphasising the see-saw mechanism, which is 
probably the most elegant way of accounting for small neutrino masses. We have further 
shown that the see-saw mechanism may be successfully applied to the atmospheric and 
solar data to yield a neutrino mass hierarchy and two large mixing angles in a technically 
natural and elegant way using the idea of sequential dominance. Sequential dominance 
requires certain mild conditions to apply, and we have seen that these conditions may in 
turn arise from the symmetries of realistic models. 
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The problem of neutrino masses and mixings should be addressed in the wider context 
of the problem of all quark and lepton masses and mixing angles, and in this wider context 
we have emphasised ideas such as unification and family symmetry which will surely play 
a role in the ultimate solution to the problem of flavour. In Table ^ we have classified 
successful models according to the different unification and family symmetries upon which 
they are based. It remains to be seen if any of these models will turn out to provide the 
solution to the problem of flavour. If this turns out to be not the case, then the effort 
will not have been in vain, since it is quite likely that some of the ideas on which these 
models are based will survive. Here we have emphasised particularly promising ideas such 
as the see-saw mechanism, sequential dominance, super symmetry, unification and family 
symmetry, which when combined with the neutrino data could help to unlock the whole 
mystery of flavour. 
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A. Equivalence of different parametrisations 

In this appendix we exhibit the equivalence of different parametrisations of the lepton 
mixing matrix. A 3 x 3 unitary matrix may be parametrised by 3 angles and 6 phases. We 
shall find it convenient to parametrise a unitary matrix by ^: 

= P2R23R13P1R12P3 (A.l) 

where Rij are a sequence of real rotations corresponding to the Euler angles 6ij, and Pi are 
diagonal phase matrices. The Euler matrices are given by 

/l 



where c 



cos 6ij and Sij 



R23 — C23 S23 
\0 -S23 C23 

Ci3 Si3 

Ri3= 10 

\ -Sl3 Ci3 

C12 S12 

R12 = -Sl2 C12 

V 1 
sinfty. The phase matrices are given by 

/l 0^ 



(A.2) 



(A.3) 



(A.4) 



P 



e'^ 
1 



(A.5) 



^It is convenient to define the parametrisation of V'' rather than V because the lepton mixing matrix 
involves V"^ ^ and the neutrino mixing angles will play a central role. 
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e*'^2 
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(A.6) 



(A.7) 



Y 

By commuting the phase matrices to the left, it is not difficult to show that the 



parametrisation in Eq.A.l is equivalent to 

yt = PU23UM2 



where P = P1P2P3 and 
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C23 S23e-'^''-' 






-S23e*'^23 
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CIS si3e-*^i3 


Ui3 = 
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5136^^" ci3 
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ci2 si2e-'^'^ 


Ui2 = 




si2e'^'^ ci2 
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1 



where 



(^23 = X + ^2 - ^3 
<Jl3 =UJi-UJ3 
612 = OJi- UJ2 



(A.8) 
(A.9) 
(A.IO) 

(A.ll) 

(A.12) 
(A.13) 
(A.14) 



The matrix U is an example of a unitary matrix, and as such it may be parametrised by 
either of the equivalent forms in Eqs.AJ or A.8 . If we use the form in Eg. A. 8 then the 
phase matrix P on the left may always be removed by an additional charged lepton phase 
rotation /SV^'^ = P\ which is always possible since right-handed charged lepton phase 
rotations can always make the charged lepton masses real. Therefore U can always be 
parametrised by 

U = U23U13U12 (A.15) 

which involves just three irremoveable physical phases 5ij. In this parametrisation the 
Dirac phase 5 which enters the CP odd part of neutrino oscillation probabilities is given 

by 

(5 = 5^3 - ^23 - 5l2- (A. 16) 



Another common parametrisation of the lepton mixing matrix is 

U = R23U13R12P0 



(A.17) 
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where 



/ e^^i \ 
e*^2 

1 



(A.18) 



and in Eq.A.17 U13 is of the form in Eg, A.10| but with 613 replaced by the Dirac phase 5. 
The parametrisation in Eq.A.17 can be transformed into the parametrisation in Eq. A.15 



by commuting the phase matrix Pq in Eq, A.17 to the left, and then removing the phases on 
the left-hand side by charged lepton phase rotations. The two parametrisations are then 
related by the phase relations 

623 = P2 (A. 19) 



J12 



/3i-/32 



(A.20) 
(A.21) 



The use of the parametrisation in Eq.A.17 is widespread in the literature, however for the 



reasons discussed in the next sub-section we prefer to use the parametrisation in Eq. [A.15 
which is trivially related to Eq. A.17| by the above phase relations. 



B. Three Family Oscillation Formulae 

At a Neutrino Factory it is relatively straightforward to measure the angle ^13 using the 
Golden Signature of "wrong sign" muons. The effect relies on the full three family oscilla- 
tion formulae which we discuss in this Appendix. For example suppose there are positive 
muons circulating in the storage ring, then these decay as /i"*" — > e'^Vf.Vfj^ giving a mixed 
beam of electron neutrinos and muon anti-neutrinos. The muon anti-neutrinos will interact 
in the far detector to produce positive muons. Any "wrong sign" negative muons which 
may be observed can only arise from the neutrino oscillation of electron neutrinos into 
muon neutrinos with probability given by a CP conserving part and a a CP violating 
part P~ . The exact formulae in vacuum are given by: 

P{ye ^ I^m) = P^i^e ^ yf.) + P'il^e ^ (B.l) 
P{i?e ^ = P^{l>e ^ + P"(l>e ^ V^) (B.2) 

where the CP conserving parts are 

P+{ue ^ i^m) = P^i^e -^^^) = - 4Re{UeiU;^U:2U^2)sm\l.27Aml^L/E) 

- 4Re{UeiU;,U:^U^3) sm\l.27AmlL/E) 

- 4Re{Ue2U;2U:3U^3)sm\l.27Aml^L/E) (B.3) 

and the CP violating parts are 

P^ {Uf, u^) = —P^{i>(, u^) = —Ci3 sin 2^13 sin 2(?i2 sin 2^23 sin 
X sin(1.27Am^iL/S)sin(1.27Am|iL/^)sin(1.27Am|2V^) (B.4) 
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Note that P~ requires all three families to contribute, and it vanishes if any mixing angle 
or mass splitting is zero. The angle ^13 may easily be extracted from Ues in the dominant 
CP conserving term P"^. 

In order to determine the CP violating phase sin 5 it is necessary to measure the CP 
violating term P~ . In order to do this one must compare the result for P(z^e ^^i) 
to the result to the case where the positive muons in the storage ring are replaced by 
negative muons and the analagous experiment is performed to measure P{i>e — > i^^i)- The 
CP violating asymmetry due to the CP violating phase 5 is given by 

from which we obtain 

5 P-(z/e^f/,) sin26li2sin(5 . . 2 ^/^^^ m a\ 

A = — ^i- ^ — sm(1.27Am2iL/S (B.6 

P+ive^Up) smells 

It is clear that in order to measure the CP asymmetry we require large 612 and large Am^i 
and this corresponds to the LMA MSW solution. In addition we require large sin 5. Also 
it would seem that having small ^13 enhances the CP asymmetry, however it should be 



remembered that the CP asymmetric rate P~ in Eq.B^ is proportional to sin 2^13, and so 
013 should not be too small otherwise the number of events will be too small. 

Unfortunately life is not quite as simple as the above discussion portrays. The Earth is 
made from matter and not anti-matter and so CP will be violated by matter effects as the 
neutrino beam passes through the Earth from the muon storage ring to the far detector. 
For example the matter effects will modify the formulas for P(fe — > f/^) involving ^13 and 

sin 2^13 



Am^i as follows 



sin 2^13 



2 

2- - cos 26*13 + sin^ 26*13 



where 



^ = 7.6 X 10'^ pE (B.8) 

where p is the density of the Earth in gcm^^ and E is the beam energy in GeV. The point 
is that for P{i'e v^) the sign of A is reversed. From one point of view this is good news, 
since unlike the vacuum oscillation formulae, Am^i enters linearly, not quadratically, and 
so matter effects enable the sign of the mass squared splitting to be determined in a rather 
straightforward way. 

However from the point of view of measuring sin 5 it leads to complications since the 



asymmetry in the rate in Eq. |B.5| can get contributions from both intrinsic CP violation 
and from matter induced CP violation, and the measured asymmetry is a sum of the two 
effects 

A^P = A^+ A'^'^tter (g_9) 
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Since both effects are by tfiemselves ratfier small, it will be a very difficult job to disentangle 
them, and the optimal strategy continues to be studied The optimal place to sit in 
order to observe CP violation seems to be at the peak of sin(1.27Am32-L/£^) in order to 
maximise P~ according to Eq. [B.4| (certainly we should avoid being at its node otherwise 
CP violation vanishes). In order to do this efficiently it may be desirable to have energy- 
tunable beams, and it is certainly necessary to have a good understanding of the density 
profile of the Earth. Assuming the LMA solution, the prospects for measuring CP violation 
at a Neutrino Factory are good. 



C. Charged lepton contributions to the lepton mixing matrix 



In this appendix we discuss the contribution of the charged lepton mixing angles to the 
lepton mixing matrix. The lepton mixing matrix is constructed in Eq j2.17| as a product of a 
unitary matrix from the charged lepton sector V^'^ and a unitary matrix from the neutrino 
sector V^^^ . Each of these unitary matrices may be parametrised by the parametrisation 
of V"^ in Eq.A.l. Thus we write 



]/El'^ _ tdEl tjEl tjEl pEL tjEl oEl 
V — ^2 ^23 -"-13 ^1 ^12 ^3 



(C.l) 
(C.2) 



where the Euler angles and phases are defined as in Eqs jA.2[jA.7| but now there are inde- 
pendent angles and phases for the left-handed neutrino and charged lepton sectors distin- 



guished by the superscripts ul and E^. The lepton mixing matrix from Eqs. 2. 17 , C.l, C.2 
is then 

U = P^^^ R^i^ P^^^ R^i^ R^i^ P^^^ P!^^Rlf^R\i,P'[^R'i'iP^^^ (C.3) 



As before we commute all the phase matrices to the left, then choose Pg 
the phases on the left-hand side, to leave just 



to cancel all 



(C.4) 



with independent phases and angles for the left-handed neutrino and charged lepton sectors. 



in the convention of Eqs.A.9, A.IO, A.ll. The phases in Eq.C.4 are given in terms of the 



phases in Eqs. C.l, C.2 by 
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(C.5) 
(C.6) 
(C.7) 
(C.8) 
(C.9) 
(C.IO) 



The form of U in Eq.|C.4| is similar to the parametrisation in Eq. 2.25 , which is the practical 



reason why we prefer that form rather than that in Eq.A.17 
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We now discuss the lepton mixing matrix to leading order in ^13 . From Eqs . A. 15 , A.9| , A.IC , A.ll , 
we find to leading order in ^13 that U may be expanded as: 



U : 
( 



C12 
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-i<5i 



^136 



-i<5i. 



-Sl2C23e^^ - Ci2S23t7l3e 
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i(<5i3— <5i2) 



C23 / 



(C.ll) 



For = 0.1, close to the CHOOZ limit, the approximate form in Eg. C.ll is accurate to 
1%. 

We now wish to expand the MNS matrix in terms of neutrino and charged lepton 



mixing angles and phases to leading order in small angles, using Eq. C.4. In technically 
natural theories, based on right-handed neutrino dominance, the contribution to ^23 comes 
mainly from the neutrino sector, ^23 ~ ^23- Furthermore in natural theories we expect 
that the contributions to ^13 are all separately small so that the smallness of this angle 
does not rely on accidental cancellations. Clearly this implies that 9^i^ and 6^^ must both 
be < ^13. Since the 13 element of U also receives a contribution from the charged lepton 
sector proportional to s^^Sgf, the same argument also implies that ~ ^13 • Therefore 
the natural expectation is that all the charged lepton mixing angles are small! Expanding 
EqjOl to leading order in small angles ^f/, 6^^ , 9^-, , we find 
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(C.12) 



where we have dropped terms of order 62^61^. The first matrix on the right hand side 
of Eq.|aT^ gives the contribution to the lepton mixing matrix from the neutrino mixing 



angles and phases, and is of the same form as Eq. C.ll . The subsequent matrices give the 
corrections to the lepton mixing matrix from the charged lepton mixing angles ^ ^13 ^ 
and efa^. 
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D. Analytic Approach to Diagonalising Mass Matrices 



D.l Proceedure for diagonalising hierarchical mass matrices 

In this appendix we discuss the diagonalisation of a general complex hierarchical matrix 
m, assuming two large mixing angles and one small mixing angle, to leading order in the 
small mixing angle, where 

^ mil mi2 mi3^ 
m = m2i m22 "1-23 
\m3i m32 m33 J 

The matrix m is diagonalised by a sequence of tranformations: 



Pi 



L ^ tdL* pL ^ pL ^ JdL*^ pR TjR pR pR pR pR 
^13 -"-23 ^2 "1^2 ^23^13M -"-12^3 



Rl2 Pi 



f mi 

7712 
^ 7713 



(D.2) 



In the case of the charged lepton mass matrix, all the rotation angles are small, while in 
the case of the neutrino mass matrix it is symmetric. The results for the general complex 
matrix m will be sufficiently general to allow us to apply them to both of the physical cases 
of interest as limiting cases. 

The proceedure for diagonalising a general hierarchical matrix m involves the following 
steps. 

1. The first step involves multiplying the mass matrix m by the inner phase matrices 



P2 defined in Eq.A.6 



Pn^ mPr 
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m2ie~*' 
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(D.3) 



m^Y m 



32 



m 



33 



The purpose of this re-phasing is to facilitate steps 2,3 using real rotation angles ^23; ^137 
as we shall see. 



2. The second step is to perform the real rotations R23 defined in Eq.A.2 on the 



re-phased matrix from stepl. The purpose is to put zeroes in the 23,32 elements of the 
resulting matrix: 



R 



23 



70,11 "^12 "^13 
™-21 "1-22 ^23 



"^31 "^32 



777 



33 



pR 



^ mil rhl2 77113 ^ 
77721 "722 

Y 77731 7773 J 



(D.4) 



The zeroes in the 23,32 positions are achieved by diagonalising the lower 23 block, using 
the reduced matrix i?23 obtained by striking out the row and column in which the unit 
element appears, to leave a 2 x 2 rotation. 
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which imphes 
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The requirement that the angles ^23 are real means that the numerators and 

denominators must have equal phases, and this is achieved by adjusting the relative phases 
(pf' — (f)^ which appear in the lower block of Eq. D.3| . The remaining elements are then given 
by the reduced rotations 

fhi2 77113 ) = ( "7l2 "^13 ) -^23 (^-^ 




R. 
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m. 



23 



21 



777, 
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31 



3. The third step is to perform the real small angle rotations i?i3 defined in Eq.A.3 on 
the matrix from step 2. The purpose is to put zeroes in the 13,31 elements of the resulting 
matrix: 

/ 77711 ^-12 

(D.IO) 



-^13 



77711 "1-12 "1-13 
77721 "^22 

y 77731 7773 



^13 



77721 "1-22 
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The zeroes in the 13,31 positions are achieved by diagonalising the outer 13 block, using 
the reduced matrix Ri^ obtained by striking out the row and column in which the unit 
element appears, to leave a 2 x 2 rotation. 



T I mil "713 
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which implies 



(D.ll) 



(D.12) 
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The requirement that the angles 9i^ and 9^^ are real fixes the absolute value of the phases 
4>^' + 4>j', since only the relative phases were fixed previously. This uses up all the phase 
freedom and thus all the resulting mass matrix elements in Eq. p.lC remain complex. Note 
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that Eq. D.IO is written to leading order in the small angles ^13, and as discussed previously 
the 23,32 elements remain zero to this order. The large complex element 7773 is approx- 
imately unchanged to this order. Due to the zeroes in the 23,32 position of the matrix 
the elements 77712 and 777-21 are also unchanged to leading order. The element 777-22 is also 
unchanged of course since it is not present in the reduced matrix. The only new element 
is therefore 

7771377731 



77711 ~ mil 



rUn 



(D.14) 
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4. The fourth step involves multiplying the mass matrix resulting from Eq. D.IO by the 

Cq 



phase matrices Pi defined in Eq.A.5 
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The purpose of this re-phasing is to facilitate step 5 using real rotation angle 6*12. 



5. The fifth step is to perform the real rotations -R12 defined in Eq. [A.4| on the re-phased 
matrix from step 4. The purpose is to put zeroes in the 12,21 elements of the resulting 
matrix: 
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The zeroes in the 12,21 positions are achieved by diagonalising the upper 12 block, using 
the reduced matrix R12 obtained by striking out the row and column in which the unit 
element appears, to leave a 2 x 2 rotation, 
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which implies 
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The requirement that the angles 6^2 and 0^ are real means that the numerators and 



denominators must have equal phases, and this is achieved by adjusting the phases xl, XR- 
6. The sixth step involves multiplying the complex diagonal mass matrix resulting 
from Eq. D.16| by the phase matrices P3 defined in Eq.A.7: 
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The result of this re-phasing is a diagonal matrix with real eigenvalues. In the case of 
charged leptons this last step can be achieved by a suitable for any choice of P^. This 
freedom in P^ enables three phases to be removed from the lepton mixing matrix. 

D.2 Diagonalising the hierarchical neutrino mass matrix 



In this appendix we shall apply the results of appendix [D.l| to the case of the complex 
symmetric hierarchical neutrino mass matrix of the leading order form of Type lA as 
shown in Table 1, which will be written in full generality as 
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where it should be remembered that for a Type lA matrix the elements in the lower 23 
block are larger than the other elements. 



The proceedure outlined in appendix D.2 for diagonalising m'^]^ is to work our way 
from the inner transformations to the outer transformations as follows. 



1. Re-phase using the • 

2. Put zeroes in the 23=32 positions using -R23- 

3. Put zeroes in the 13=31 positions using Ri^. 

4. Re-phase the mass matrix using P^^ . 

5. Put zeroes in the 12=21 positions using -R^f- 



Make the diagonal elements real using the P^^ . 



1 "^1 



However in the hierarchical case rruk , mXk are smaller than m,^ 



If is small, then for the hierarchical case ^ m2 this proceedure will result in an 
approximately diagonal matrix to leading order in . One might object that after step 
3 the i?^3 rotations will "fill-in" the zeroes in the 23,32 positions with terms of order 
multiplied by 7x1^2 , 
by a factor of 6'^^ 

factor of (^if)^ compared to the 33 element. This means that after the 5 steps above a 
hierarchical matrix will be diagonal to leading order in 9'^^, as claimed. For the inverted 
hierarchical neutrino case a different proceedure must be followed, as discussed in the next 
sub-section. Here we shall systematically diagonalise the hierarchical neutrino mass matrix 



which means that the "filled-in" 23,32 entries are suppressed by a total 



in Eq.D.21 by following the above proceedure as follows. 



The first step is to re-phase the matrix in Eq.D.21 using P^^* so that the neutrino 
mass matrix becomes. 
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To determine the 23 neutrino mixing angle 



^23 ^® perform a 23 rotation which diagonalises 



the lower 23 block of Eq. p.22| . From Eq. p.q we find the 23 neutrino mixing angle ^23 



The relative phase 
real. 
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is fixed by the requirement that the angle ^23 Eq.D.23 be 
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(DM) 



After the 23 rotation in Eq.D.4, the neutrino mass matrix in Eq.D.22 becomes 
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The lower block elements are given by 
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which implies 
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and from Eq.l 
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We now perform a 13 rotation on the neutrino matrix in Eq. p.25| which diagonalises 

13 



the outer 13 block of Eq.D.25 and determines the 13 neutrino mixing angle O'^k. From 



Eq jD.12| we find the 13 neutrino mixing angle as 

^13 



^13 



(D.30) 



The absolute phases 02^, 4>-^^ are fixed by the requirement that the angle O-^f^ in Eq. D.30 
be real, 



S23 I"^i2l sm((Ai2 - 02"" - (t>'z) + C23 I"^i3l sin 
After the 13 rotation in Eq. D.lC , Eq. D.25| becomes 
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To leading order in the only new element in Eq. D.32| is 
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It only remains to determine the 12 neutrino mixing angle by diagonalising the 
upper 12 block of Eq p32[ From Eq. plSj we find the 12 neutrino mixing angle 6*^2 



tan 29^^ 
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The phase x*^^ is fixed by the requirement that the angle ^J'o in Eq.D.34 be real, 



\m22 \ sm' 



22 



<^r2-x''") = imriisin(0ri-0r2+x'^") 



(D.35) 



-51 - 



After the 12 rotation the upper block of the matrix in Eq. D.32 is diagonal and the resulting 
matrix is 

/ m'l \ / mie'^i \ 

m2e'^2 (D.36) 
mse^'^a / 
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m'2 I = 
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where from Eq.D.17 
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It is a simple matter to adjust the phases in P^^ to remove the phases in Eg. p. 36 



and make the neutrino masses real, as in Eq.D.2C, 



(D.39) 



This completes the diagonalisation. In the case of neutrino masses, unlike the case of the 
charged fermions, there is no left over phase freedom. This is the reason why the lepton 
mixing matrix has three more physical phases than the CKM matrix. 

D.3 Diagonalising the inverted hierarchical neutrino mass matrix 

In this appendix we shall consider the case of the complex symmetric inverted hierarchical 
neutrino mass matrix of the leading order form of Type IB in Table 1. In this case the 
proceedure is as follows. 

1. Re-phase m^^;^ using the ■ 

2. Put zeroes in the 13=31 positions using B!^^. 

3. Put zeroes in the 23=32 positions using R^k^. 

4. Re-phase the mass matrix using P^^. 

5. Put zeroes in the 12=21 positions using R\2- 

6. Make the diagonal elements real using the P3 ^ . 



We continue to write the neutrino mass matrix as in Eq jD.2l| , but now it should be 
remembered that for a Type IB matrix the 12,13 elements are now larger than the other 
elements. This is the reason why the above proceedure differs from that for the case of the 
hierarchical neutrino mass matrix. 
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We first perform the re-phasing as in Eq. D.22 . Then we determine the 23 neutrino 
mixing angle by performing a 23 rotation such that 
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which gives the 23 neutrino mixing angle ^23 ™ this case to be 
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Since the Euler angles are constrained to satisfy 
and this then fixes 
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The lower block elements are given by 
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and the remaining diagonal elements are given as before in Eqs. p^p!28| 
We next perform a small angle 13 rotation such that 
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Note that to leading order in 0^3 the large element is unchanged. The remaining 

h3 ■ 



elements in Eq. D.48 are also unchanged to leading order in 6'^!^. The 13=31 element in 
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Eq, D.48 gets filled in by a term 0'^i^{m' 



V 

12 



7713) which is of order (^^f )^ compared to m 
and does not appear to leading order in Q^. From Eq. D.4S| the 13 neutrino mixing angle 
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The requirement that is real fixes the absolute value of the phases (^2^ , 



The left hand side of Eq.D.48 now resembles the left hand side of Eq.D.32, except that 
here is unchanged due to the zero 13=31 element after the 23 rotation. Therefore 
the rest of the diagonalisation process follows that of the previous hierarchical case from 



Eq,D.34 onwards, where now 



tan 29^^ 



Note that in the inverted hierarchy case here we have 



2 






1 771^2 







(D.51) 
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which implies an almost degenerate pair of pseudo-Dirac neutrino masses (with opposite 
sign eigenvalues) and an almost maximal 12 mixing angle from Eq.D.51. 
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